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NOMENCLATURE 


a - Superpressure  at  nadir  of  balloon  (feet) 

2 

A - Cross-sectional  area  in  horizontal  plane  (feet) 

b - Specific  lift  defined  as  difference  in  specific  weights  of  atmosphere 
and  inflation  gas  (pounds/feet  ) 

d(  ) - Differential  quantity  as  defined 

D - Constant  used  to  characterize  film  behavior  as  defined  by  equations 
[12]  and  [13]  (Stress)"' 

F - Tensile  force  carried  by  an  individual  load  tape  (pounds) 

Ky  - Constant  used  to  characterize  tape  behavior  as  defined  by  equation  [14] 

(pounds) 

l - Gore  width  at  any  gore  position  (feet) 

N - Number  of  gores 

o 

p - Pressure  differential  across  film  (pound/feet  ) 

r - Radial  location  of  load  tape  relative  to  balloon  centerline  in  a 
horizontal  plane  (feet) 

R - Radius  of  curvature  of  film  or  tape  (feet) 

s - Gore  position  measured  from  nadir  (feet) 

T - Total  force  in  the  meridional  direction  divided  by  2n  (pounds) 

V - Volume  of  the  balloon  (feet^) 

2 

W - Weight  per  unit  area  or  length  of  film  or  tape  respectively  (Ibs/ft  or  Ibs/ft) 

z - Vertical  height  of  an  arbitrary  point  measured  from  the  nadir  (feet) 

a - Angle  in  plane  containing  circumferential  radius  of  curvature  defined 
by  equation  [28] 

e - Angle  in  plane  containing  circumferential  radius  of  curvature  shown  in 
Figure  3. 

1 - Angle  defining  circumferential  radius  of  curvature  shown  in  Figure  3. 

A(  ) - Symbol  denoting  a change  in  quantity 

e - Normal  strain  defined  as  the  change  in  length  per  unit  length 

e - Angle  made  by  the  tangent  to  the  load  tape  in  the  meridional  direction 
and  the  balloon  centerline. 

2 

a - Normal  stress  or  stress  resultant  as  defined  (pounds/feet  ) 

2 

t - Shear  stress  or  stress  resultant  as  defined  (pounds/feet  ) 

SUBSCRIPTS 

c - Circumferential  direction 
d - Design  Value 
f - Film  value 

i - Arbi  bra ry  point  in  finite  difference  equation 
m - Meridional  direction 
o - Undeformed  value 
T - Tape  value  7 


I.  INTRODUCTION 


The  design  and  analysis  of  free  balloons  for  scientific  application  has 
received  the  attention  of  many  investigators  over  the  years.  The  pioneering 
efforts  by  Upson  (1)  and  the  University  of  Minnesota  (2)  contributed  to  the  well 
documented,  computerized,  shape  and  stress  calculations  of  Smalley  (3,4)  in  both 
the  fully  deployed  and  partially  deployed  state.  This  effort,  although  completed 
over  a decade  ago,  continues  to  be  the  foundation  on  which  successful  balloon  de- 
signs are  generated. 

The  investigations  by  manufacturers  and  research  organizations  in  recent  years 
have  been  devoted  to  the  application  of  these  procedures  to  produce  various  designs, 
attempts  to  characterize  the  film  material,  or  detailed  analysis  of  the  state  of 
stress  and  deviations  from  the  design  shape  due  to  material  deformation.  The  most 
sophisticated  of  the  analysis  techniques  utilizes  finite  elements  and  the  minimi- 
zation of  potential  energy  to  obtain  deviations  from  the  design  shape.  Due  to  the 
nonlinearities  associated  with  large  displacements,  the  adaptation  of  such  a pro- 
gram (5)  to  include  the  loads  associated  with  inflatables  is  a significant  achieve- 
ment. The  effects  of  lobing  have  been  analyzed  by  Alexander  (6)  in  a manner  which 
incorporates  the  characterization  of  polyethelene  films  with  the  design  shape.  In 
addition,  Rand  (7)  has  utilized  finite  difference  techniques  to  determine  the 
effects  of  load  introduction  into  pressurized  films. 

Many  attempts  have  been  made  to  characterize  candidate  films  and  develop  the 
necessary  testing  apparatus  to  evaluate  their  material  properties.  Considerable 
effort  has  been  expended  in  an  attempt  to  evaluate  the  cold  brittleness  of  a film 
as  a measure  of  the  quality  of  the  material.  Both  uniaxial  and  biaxial  results 
have  been  reported  by  Weissmann  (8),  Alexander  (9),  and  Kawada  (10),  et  al.  In 


fr 


addition,  the  characterization  of  reinforced  films  has  been  studied  by  Alley  (11) 
and  Munson  (12).  Their  results  indicate  that  in  the  future,  composite  films  may 
be  manufactured  with  overall  properties  which  are  optimum  for  a particular  design. 
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I 
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’! 
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Each  of  the  studies  previously  mentioned  as  well  as  manufacturing 
improvement  in  film  and  seal  quality  have  been  sound  engineering  programs. 

As  a result  many  assumptions  have  been  made  regarding  the  shape,  stress  and  failure 
mechanisms  operative  in  balloon  systems.  These  assumptions  have  been  necessary  to 
obtain  operational  systems  when  they  were  needed.  However,  Dwyer  (13)  has  attempted 
to  quantify  the  capability  for  growth  of  free  balloons  in  a manner  which  suggests 
that  the  state-of-the-art  may  have  matured  to  the  point  of  limiting  growth.  This 
indicates  that  there  is  a need  to  reassess  those  assumptions  that  were  made  to 
achieve  the  current  state-of-the-art  in  the  hope  of  developing  an  innovative  design 
concept. 

At  the  present  time  the  design  process  and  stress  analysis  are  carried  out  in- 
dependently. This  is  necessitated  by  the  fact  that  the  shape  must  be  known  apriori 
in  order  to  perform  any  of  the  analyses  previously  mentioned.  The  finite  element 
analysis  technique  requires  the  coordinates  of  each  node  of  every  element  prior 
to  initiation  of  the  solution  scheme.  Although  very  detailed  solutions  can  be 
obtained  by  this  technique,  in  order  to  interface  such  a technique  with  an  existing 
design  program  would  require  an  iterative  technique  to  be  developed  and  excessive 
amounts  of  computer  time.  The  analysis  technique  of  Alexander  (6)  assumes  that 
the  shape  predicted  by  the  design  program  is  in  fact  the  shape  of  the  gore  seams. 
This  is  somewhat  presumptuous  since  any  limitation  of  the  design  shape  is  auto- 
matically included  in  the  stress  analysis. 

The  analysis  of  Rand  (7)  in  determining  the  stresses  associated  with  load 
introduction  into  a pressurized  film  involves  the  simultaneous  solution  of  two 
partial  differential  equations  in  terms  of  two  unknown  displacements  at  each  of 
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200  points.  The  technique  suffers  from  the  same  limitations  as  the  finite  element 
routine  in  that  it  would  require  an  initial  shape,  an  iterative  technique  and  ex- 
cessive computer  time. 


The  initial  shape  calculations  on  which  all  analyses  are  based  are  not  without 
certain  simplifying  assumptions  which  make  the  problem  tractable.  By  assuming  that 
the  meridional  stress  is  uniformly  distributed  around  a symmetrically  deployed  mem- 
brane, the  two  equations  of  equilibrium  become  sufficient  for  a unique  solution. 
Unfortunately,  the  stresses  are  not  uniformly  distributed  throughout  the  membrane 

i 

but  instead  are  transferred  from  the  load  tapes  to  the  gore  panels  through  a vari- 
ety of  mechanisms.  It  is  presumed  that  near  the  apex  and  nadir  of  a free  balloon 
in  the  fully  deployed  configuration,  the  meridional  loads  are  carred  almost  en- 
tirely by  the  load  tapes  while  near  the  maximum  radius  of  the  balloon  this  load 
is  carried  almost  entirely  by  the  film  A similar  situation  will  exist  during  the 
ascent  phase  of  the  flight  except  that  excess  material  must  be  carried  as  pay- 
load  and  the  mechanism  of  load  transfer  may  be  altered. 

It  has  been  noted  by  Smalley  (14)  and  others  that  the  capability  of  placing 
a specified  payload  at  a specified  altitude  is  limited  only  by  the  strength  to 
weight  ratio  of  available  films  and  the  ability  to  seal  and  handle  the  film.  It 
should  be  noted  that  as  efforts  continue  by  the  manufacturers  to  improve  the 
material  and  handling  properties  of  available  films  and  candidate  films,  the 
development  of  a higher  order  design  procedure  may  lead  to  lighter  vehicles  with 
the  same  capabilities  as  the  present  generation  of  vehicles. 


II.  DESIGN  FORMULATION 


The  normally  accepted  equations  and  assumptions  for  the  shape  and  stresses 
in  a balloon  have  been  thoroughly  documented  by  Smalley  (3).  The  assumptions 
are  as  follows: 

1.  The  balloon  is  rotationally  symmetric  about  a vertical  axis. 

2.  Meridional  and  circumferential  stresses  are  constant  at  all  points  on 

the  circle  formed  by  an  intersecting  plane  normal  to  the  axis  of  symmetry. 

This  preculdes  the  possibility  of  shear  in  the  balloon. 

3.  The  densities  of  the  inflation  gas  and  surrounding  air  are  constant. 

4.  The  balloon  material  is  inextensible,  thin,  and  incapable  of  supporting 
any  bending  or  compressive  loads. 

The  resulting  set  of  differential  equations  include  two  equilibrium  equations 
and  four  geometric  relationships  which  must  be  solved  simultaneously  for  the 
shape  and  stress  variables  at  each  point  along  the  gore.  This  approach  is  known 
in  mechanics  as  the  "stress  formulation"  which  is  the  preferable  method  when  all 
of  the  boundary  conditions  may  be  expressed  in  terms  of  forces  rather  than  dis- 
placements. The  resulting  shape  from  this  formulation  is  the  deformed  configuration 
but,  since  the  film  is  assumed  to  be  inextensible  then  the  undeformed  shape  and  the 
manufactured  shape  will  be  identical. 

The  meridional  stress  is  computed  as  one  of  the  variables  while  the  circum- 
ferential stress  must  be  known  in  at  least  some  functional  form  prior  to  attempting 
a solution.  The  most  common  shape  is  know  as  the  "natural  shape"  which  is  the 
configuration  that  results  when  th°  circumferential  stress  is  zero.  Other  designs 
may  be  obtained  by  allowing  the  circumferential  stress  to  be  proportional  to  the 
meridional  stress,  or  the  change  in  meridional  stress  or  any  other  functional  relation- 
ship. This  is  an  acceptable  technique  since  different  shapes  will  result  from 
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the  different  assumed  stress  relationships.  The  most  well  know  configurations 
are  those  natural  shapes  listed  in  the  Sigma  tables  by  Smalley  (3).  The  stresses 
predicted  by  this  method  would  be  accurate  If  the  balloon  could  be  built  in  a 
manner  which  obeys  the  various  assumptions  made  in  the  formulation  of  the  equations. 
Unfortunately,  it  would  be  necessary  to  build  the  complete  balloon  as  a shell  of 
revolution  without  seams  or  load  tapes  in  order  to  achieve  this  shape  and  stress. 

Another  disconcerting  result  of  these  assumptions  is  that  the  meridional  film 
load  is  formulated  as  the  product  of  the  radius  (measured  in  the  horizontal  plane) 
and  the  meridional  stress.  This  product  remains  finite  throughout  the  balloon 
but  as  the  radius  approaches  zero  near  the  ends  of  the  balloon,  the  meridional 
stress  necessarily  approaches  infinity.  This  value  of  stress  is  usually  ignored 
as  a singularity  associated  with  the  coordinate  system  and  the  practical  realization 
that  the  load  tapes  will  prevent  such  stresses  from  developing.  Another  mitigating 
factor  is  the  fact  that  film  thickness  is  not  dictated  by  the  stress  predicted 
by  this  procedure  but  by  practical  manufacturing  limitations.  However,  as 
balloons  become  larger  and  payloads  become  heavier,  the  normally  safe  designs 
may  become  marginal. 

. 

In  order  to  obtain  a more  realistic  estimate  of  the  state  of  stress  the 
design  equations  have  been  rederived  in  a manner  to  take  into  account  the 
effects  of  load  tape  stiffness  and  film  modulus.  While  the  equations  used  by 
Smalley  were  derived  for  a differential  element,  the  following  formulation  must 
be  considered  in  an  integral  sense.  It  is  again  assumed  that: 

1.  The  balloon  is  rotationally  symmetric  about  a vertical  axis. 

2.  Meridional  and  circumferential  stresses  are  constant  at  all  points  on 

the  circle  formed  by  an  intersecting  plane  normal  to  the  axis  of  symmetry. 

3.  The  densities  of  the  inflation  gas  and  surrounding  air  are  constant. 

4.  Both  film  and  tape  material  are  linearly  elastic  and  orthotropic. 
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Balloon  Coordinates 


F+dF 


(b)  Forces  Acting  on  Di crete  Element 


Figure  1 - Discrete  Balloon  Design  Element 


5.  The  meridional  strain  in  the  film  is  equal  to  the  meridional  strain  in 


the  tape. 

The  discrete  element  is  shown  in  Figure  1 and  represents  a portion  of  the 
surface  bounded  by  two  horizontal  planes  separated  by  a distance  dz  and  two 
vertical  planes  containing  the  axis  of  symmetry  and  the  centerlines  of  two 
adjacent  gores.  The  angle  separating  these  two  planes  is  then  2*/N  where  N is 
the  number  of  gores.  By  summing  forces  in  the  meridional  direction  it  can  be 
shown  that: 

^d(ram)  + dF  = sineds  + rW^  + W^Jcoseds 

d(r°m  + 77)  nWt 

Therefore,  g-$--  - ocsine  + (rWf  + -jocose 


[1] 


After  summing  forces  perpendicular  to  the  element  the  following  equation  results: 
NFx  de  _ . NW 


<% + sr>  3?  ■ °ccose  ’ <r“f  * ^,s1ne ' pr 

If  the  meridional  load  parameter  is  defined  as: 

T - „ + NF 

T = ro_  + 9— 

m ZTT 

And  if  the  distributed  weight  parameter  is  defined  as: 

rW  = rW  + — ^-T 
rw  _ rwf  + 2it 

Then  equations  [1]  and  [2]  may  be  rewritten  in  the  following  form: 


[2] 


[3] 


[4] 


JT 

= ocsine  + rWcose 
d6 

and  T^j  = accose  - rWsine  - br(z  + a) 

Where  the  pressure  differential  is  given  as  usual  by: 
p = b(z  + a) 

The  usual  geometric  relations  are: 


[5] 

[6] 

[7] 


£ ■ «i» 

dz 

d?  = cose 
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[8] 

[9] 
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dA 

d? 


[10] 


1 

s 

ii 
; ■ 

i 

; 

i 

B 


I 


= 2 nr 

gj  = nr^cose  [11] 

In  this  form  the  set  of  equations  Is  identical  to  those  used  by  Smalley; 
therefore,  any  existing  design  may  be  used  to  obtain  the  variables  in  question, 
i.e.,  T,  9,  r,  z»  A,  and  V.  However,  the  difference  between  this  procedure  and 
previous  approaches  to  the  problem  is  in  the  evaluation  of  the  state  of  stress 
and  the  manufactured  geometry  of  the  balloon. 


III.  DESIGN  PROGRAM  DEVELOPMENT 


l A 


The  essential  difference  between  the  usual  design  procedure  and  the  procedure 
described  here  lies  in  the  interpretation  of  the  final  results.  The  total  merid- 
ional load,  equation  [3],  must  be  partitioned  between  the  film  and  tapes  according 
to  their  respective  material  properties.  When  this  is  done,  the  resulting  film 
stresses  may  be  used  in  conjunction  with  the  appropriate  constitutive  equation 
for  the  film  to  determine  the  undeformed  or  stress  free  shape  of  the  balloon. 

This  is  the  shape  that  should  be  used  for  manufacturing  purposes  instead  of  the 
deformed  shape  as  is  presently  used. 

If  it  is  assumed  that  the  film  is  linearly  elastic  and  orthotropic,  then 
the  governing  constitutive  equations  may  be  written  as: 


£ = Da  + Dma_ 
m mm  me  c 


e _ D a + D a 
c = me  m c c 


[12] 

[13] 


It  should  be  recalled  that  the  circumferential  stress  is  known  as  an  input  to 
the  design  process.  As  a common  example,  the  fully  tailored  natural  shape  design 
assumes  a to  be  zero.  Equation  [3]  may  now  be  rewritten  assuming  that  the  merid- 
ional  strain  in  the  film  and  tape  are  equal. 


Here  Kj  is  the  tape  modulus  and  relates  the  tape  force  to  the  tape  strain  in  a 
linear  fashion. 

The  undeformed  shape  of  the  balloon  may  now  be  found  from  the  definition  of 
strain.  In  particular,  the  original  radius  of  the  balloon  at  any  gore  position 
is  found  from  the  circumferential  strain: 


[16] 


as  rosvrrrj 

where  ec  is  found  from  equations  [13]  and  [14].  The  undeformed  gore  position  is 
determined  from  the  definition  of  meridional  strain.  In  particular: 

Therefore 


The  original  gore  position  is  then  found  by  itegration  of  equation  [18]  to  any 
deformed  position,  s,  along  the  gore. 

In  order  to  perform  these  computations  it  is  necessary  to  have  an  efficient 
routine  design  program  that  will  accomodate  a variety  of  films,  tapes  and  flight 
conditions.  Therefore,  a balloon  design  program  was  written  on  the  basis  of  the 
"BALLOON"  and  "FLATBALL"  programs  reported  by  Smalley  (3).  These  programs  were 
modified  somewhat  to  include  an  automatic  computation  of  the  starting  angle  which 
significantly  speeds  convergence.  The  program  consists  of  a variety  of  subroutines 
which  eliminates  the  need  to  manually  compute  a variety  of  parameters  needed  for 
the  design  process.  A listing  of  this  program  is  contained  in  Appendix  A with 
sufficient  documentation  to  permit  use  of  the  program. 

The  various  subroutines  perform  the  following  functions: 

A.  BALDE-This  is  an  executive  routine  which  reads  the  required  input 
information  on  the  type  of  balloon  to  be  designed.  After  the  design 
has  been  established  it  computes  the  load-altitude  curve  of  the  balloon. 

B.  SUBROUTINE  DESIGN  - This  routine  is  based  on  the  computational  scheme 
reported  by  Smalley  (3).  It  is  modified  to  compute  the  various  stresses 
in  accordance  with  equation  [14]  and  computes  the  balloon  table  layout 


on  the  basis  of  equations  [16]  and  [18]. 
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C.  SUBROUTINE  MATL  - This  routine  computes  the  various  material  properties 
needed  by  the  design  process  such  as  film  and  tape  stiffness  (modulus) 


and  weight.  The  computations  are  based  on  data  from  a variety  of 

: 

sources  and  are  empirically  fitted  functions  of  temperature. 

D.  SUBROUTINE  ATMOS  - This  is  a program  for  computing  the  various  standard 
properties  and  specific  lift  when  the  design  altitude  is  given  in  meters. 

E.  SUBROUTINE  ATM0S.2  - This  program  performs  the  same  function  as  SUBROUTINE 
ATMOS  except  the  design  altitude  is  given  in  millibars  of  pressure. 

F.  SUBROUTINE  BOYNCY  - This  routine  uses  the  standard  atmosphere  to  compute 
the  pressure  altitude  when  the  specific  bouyancy  is  known.  It  is  used 
in  the  computation  of  the  load-altitude  curve. 

G.  SUBROUTINE  GLNGTH  - This  is  an  empirical  fit  of  the  unique  gore  length- 
gross  load  curve  reported  by  Smalley  (15). 


IV.  DESIGN  RESULTS 


In  order  to  demonstrate  the  impact  of  this  procedure  a number  of  studies 
have  been  performed  which  demonstrate  the  influence  of  various  parameters  such 
as  film  thickness  and  number  of  load  tapes  on  the  film  stress.  One  such  study 
has  been  presented  by  Keese  (16)  where  he  concluded  that: 

A.  The  total  surface  area  of  a balloon  and  its  manufactured  cost  may  be 
significantly  reduced  if  the  design  were  based  on  a film  stress  including 
load  tape  effects. 

B.  The  total  surface  area  of  a balloon  and  its  manufactured  cost  may  be 
modestly  reduced  by  simultaneously  increasing  the  number  of  load  tapes 
and  reducing  the  film  thickness. 

Although  these  conclusions  may  have  been  premature,  the  existence  of  a tolerable 
state  of  stress  is  predicted  under  design  conditions. 

As  a second  example  of  this  program,  a balloon  was  designed  which  has  a 
very  similar  set  of  characteristics  as  a 20.8  MCF  balloon  used  by  Ale*"  der  (6) 
to  study  the  effects  of  lobing.  The  balloon  was  designed  as  a fully  tailored 
natural  shape  with  the  following  design  parameters: 


Design  Altitude: 

126,000  ft. 

Design  Payload: 

2,700  lbs. 

Maximum  Payload: 

5,250  lbs. 

Film  Material : 

Polyethelene 

Shell  Thickness: 

.0008  in. 

Cap  Thickness: 

.0017  in. 

Load  Tape  Material 

: Polyester 

Load  Tape  Rating: 

400  lbs. 

Inflation  Gas: 

Helium 

Sample  input  and  output  for  this  particular  balloon  are  given  in  Appendix  B. 

The  integration  step  DSO  was  adjusted  until  the  final  increment  of  integration 
at  the  top  of  the  balloon  was  approximately  equal  to  all  other  steps.  The 
resulting  design  contained  202  points  (201  increments)  and  has  the  following 
characteristics: 

Balloon  Volume:  20.9  MCF 

Surface  Area:  375,000  ft.2 

Film  Weight:  2214  lbs. 

Tape  Weight:  561  lbs. 

Deformed  Gore  Length:  530.2  ft. 

Undeformed  Gore  Length:  527.4  ft. 

Deformed  Half  Gore  Width:  50.60  in  (maximum) 

Undeformed  Half  Gore  Width:  50.72  in  (maximum) 

It  may  be  noted  that  the  maximum  difference  between  the  deformed  and  undeformed  half 
gore  width  is  .12  inch  which  is  probably  not  within  manufacturing  tolerances. 

However,  this  small  difference  may  cause  a strain  of  0.2%  in  the  circumferential 
direction.  It  may  also  be  noted  that  this  balloon  will  elongate  2.8  feet  from 
its  manufactured  length  when  loaded  under  design  conditions.  However,  by  building 
the  balloon  according  to  its  deformed  length,  excess  material  is  included  in 
the  meridional  direction.  This  feature  may  serve  to  reduce  the  strain  in  that 
direction. 

The  meridional  stress  predicted  by  the  usual  assumption  of  a tapeless 
construction  is  shown  for  comparison  purposes  with  that  stress  predicted  by 
equation  [14]  in  Figure  2 in  non-dimensional  form.  It  may  be  noted  that  the 
film  stress  is  significantly  less  than  the  stress  normally  assumed  and  always 
remains  finite.  The  balance  of  the  forces  in  the  meridional  direction  is  carried 
by  the  load  tapes  under  these  conditions.  However,  in  both  cases  the  circumferential 
stresses  are  zero  which  is  an  inherent  assumption  in  the  natural  shape  design. 
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Meridional  Stress  Distribution 


V.  ANALYSIS  FORMULATION 


It  has  been  shown  that  the  inclusion  of  material  properties  is  relatively 
straight  forward  provided  the  design  shape  is  assumed  to  be  correct.  However, 
the  stress  analysis  of  a flexible  structure  is  significantly  different  from  the 
design  problem.  Once  the  balloon  has  been  manufactured,  the  available  material 
at  each  gore  position  is  fixed  and  the  solution  should  be  obtained  as  a boundary 
value  problem  rather  than  as  an  initial  value  problem.  In  addition,  the 
simplifying  assumptions  of  zero  circumferential  stress  and  a surface  of 
revolution  should  not  be  made  since  it  is  suspected  that  this  will  significantly 
alter  the  state  of  stress.  This  problem  has  been  recognized  for  many  years 
and  various  researchers  have  obtained  solutions  with  varying  degrees  of  success. 

Gilbert  (17)  was  able  to  successfully  model  a parachute  structure  utilizing 
a global  transformation  matrix  to  relate  the  local  coordinates  to  generalized 
coordinates.  However,  in  order  to  optimize  the  aerodynamic  drag  characteristics 
of  the  system,  he  simplified  his  model  to  the  case  of  zero  circumferential 
stress  which  would  correspond  to  the  design  shape  of  a balloon  system.  Although 
the  formulation  of  the  differential  equations  is  similar  to  the  approach  to  be 
presented  here,  the  simplifying  assumptions  and  radically  different  boundary 
conditions  prevent  the  use  of  this  program.  Alexander  solved  this  problem  for 
several  balloon  systems  by  assuming  that  the  load  tapes  assumed  the  design  shape 
while  the  film  was  permitted  to  bow  out  between  the  tapes.  He  permitted  the 
meridional  radius  of  curvature  to  change  across  the  gore  but  this  in  turn  prevented 
the  equilibrium  equation  in  the  meridional  direction  from  being  satisfied. 
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The  formulation  to  be  presented  here  will  postulate  a mechanism  by  which 
loads  are  transfered  between  the  load  tapes  and  the  film.  The  combination  of 
shear  stress  and  lobing  between  the  load  tapes  will  permit  the  coupling  of  the 
tape  shape  to  the  film  shape.  In  the  process  of  developing  this  formulation  a 
number  of  features  will  be  introduced  for  the  first  time.  In  order  to  eliminate 
as  much  confusion  as  possible,  Lagrangian  coordinates  will  be  used  whenever 
possible. 

! 

o 

Lagrangian  coordinates  are  routinely  used  in  time  dependent  problems  in  the 
area  of  fluid  and  solid  mechanics  to  denote  the  position  of  the  particles  initially. 

: ■ 

However,  this  system  of  coordinates  is  ideally  suited  to  problems  in  elasticity 
where  material  properties  have  been  defined  with  respect  to  the  undeformed 
dimensions.  All  balloon  material  properties  now  being  developed  at  Texas  A&M 
University  utilize  "engineering"  stress  and  strain  rather  than  "true"  stress 

? 

and  strain.  Therefore,  in  this  system  of  coordinates,  the  stress  is  defined 
as  the  load  per  unit  initial  (as  manufactured)  area  rather  than  the  actual  area. 

The  Lagrangian  gore  length  remains  the  manufactured  length,  regardless  of  the 
actual  deformed  length  of  the  gore.  Another  important  quantity  which  is 
simplified  by  this  technique  is  the  increment  of  balloon  film  or  tape  weight 
between  any  two  gore  positions.  This  increment  is  known  from  the  manufactured 
shape  and  remains  fixed  regardless  of  deformation. 

In  general,  it  will  be  assumed  that  each  balloon  gore  deploys  in  an  identical 
manner  about  the  balloon  centerline.  All  deformations  will  be  assumed  to  be 
linearly  elastic  and  orthotropic  although  numerical  results  have  been  obtained 
only  for  the  isotropic  case.  The  film  will  be  permitted  to  lobe  between  the 
load  tapes  resulting  in  a surface  which  is  not  a surface  of  revolution.  Due  to 
the  assumed  symmetry  of  deployment,  equilibrium  equations  will  be  written  for  a 
single  differential  element  of  tape  and  an  adjacent  element  of  film  at  an 
arbitrary  gore  position. 
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The  usual  deformed  balloon  coordinate  description  as  shown  in  Figure  3 will 
be  used  to  demonstrate  the  compatibility  of  this  formulation  with  the  familiar 
design  equations.  The  film  is  assumed  to  lobe  in  a plane  containing  the  merid- 
ional radius  of  curvature  to  two  adjacent  load  tapes.  The  deployed  distance 
of  the  load  tape  from  the  centerline  in  the  horizontal  plane  is  designated 
r. 

Consider  a differential  length  of  tape  at  some  arbitrary  point  on  the  balloon 
as  shown  in  Figure  4.  The  forces  acting  on  this  element  include  a changing  tape 
force,  F,  as  well  as  forces  due  to  the  film  stresses,  x and  aQ,  and  the  tape 
weight.  It  may  be  noted  at  this  point  that  the  actual  length  of  the  element  is 
ds  whereas  the  original  length  was  dsQ.  Since  "engineering"  stress  will  be  used 
consistently,  the  original  film  area  will  betQdsowhere  tQ  is  the  undeformed 
film  thickness.  Therefore,  the  stress  will  always  be  multiplied  by  t so  that 
the  units  of  stress  may  simply  be  considered  to  be  load  per  unit  original  length. 

Summing  forces  in  the  vertical  direction  it  can  be  shown  that: 

(F+dF)cos(e+de)  - Fcose  - Wyds0  - 2 x cose  dsQ  - 2 ac  sine  sine  dsQ  = 0 

The  circumferential  stress  term  arises  from  the  fact  that  e is  defined  in  a 
plane  containing  the  meridional  radius  of  curvature  and  not  the  horizontal  plane. 
This  equation  may  be  rewritten  by  taking  the  limit  as  dsQ  and  de  approach  zero  as: 

<Kfcos_el  _ + 2xcose  + 2acsine  sine  [19] 

o 

A second  differential  equation  may  be  obtained  by  summing  forces  in  the  hori- 
zontal plane.  The  resulting  differential  equation  after  taking  the  limit  as  above 
becomes : 

diFsiHel  - 2xsine  - 2ocsine  cose  [20] 

o 

Equations  [19]  and  [20]  express  the  equilibrium  of  forces  acting  at  any 
point  along  the  load  tape. 


1 
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Consider  now  a differential  length  of  film  that  runs  between  two  adjacent 
load  tapes  at  the  same  gore  position  as  shown  in  Figure  5.  Here  it  is  assumed 


n 


f 


that  a constant  meridional  stress  is  uniformly  distributed  across  the  gore 
width  which  is  2irrQ/N  where  N is  the  number  of  gores  and  rQ  is  the  manufactured 
radius  of  the  balloon.  The  meridional  stress  is  permitted  to  change  in  the 
meridional  direction  as  shown.  In  addition  there  will  be  forces  due  to  the 
shear  stress,  t,  the  circumferential  stress,  oc,  the  film  weight  which  is 
expressed  in  Lagrangian  coordinates  and  the  force  due  to  pressure.  This  last 
force  is  obtained  by  using  the  actual  projected  area  which  may  then  be  expressed 
in  the  original  coordinates  by  observing  the  definition  of  meridional  strain,  i.e.. 


ds 

ds„ 


= 1+ 


[17] 


Summing  forces  in  the  vertical  direction  it  can  be  shown  that: 

^(r«am+d(r«°m))cos(0+d0)  — r^-r^a  cose— + 2-rcose  dsrt 
N'  o m ' o nr  ' Norn  Nofo  o 

-2prsin(jj-)sine  ds  + 2ocsinB  sine  dsQ  = 0 

Utilizing  equation  [17]  and  taking  the  limit  as  dsQ  and  de  approach  zero,  this 
equation  becomes: 

dl^mqp.se.)  = r u + pr(l+E  )-sin(Jr)  sine  - ^cose  - \sine  sine  [21] 

ds^  of  r m it  N it  it  c 

o 

A second  differential  equation  for  the  film  may  be  obtained  by  summing  forces 
in  the  horizontal  plane  and  taking  the  limit  as  before.  The  resulting  equation  is: 

jir^ntSjP6)  = -pr(l+em)JJsin( jocose  - ^sine  + ^c(sin(J-)cose  + cosj sine  cose) 


In  working  with  these  equilibrium  equations  it  is  a simple  matter  to  rotate 
the  system  in  such  a way  that  the  equilibrium  equations  in  the  meridional  and 
normal  directions  are  obtained.  When  this  is  accomplished  the  equations  become: 
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Meridional  Tape  Equilibrium 

= 2t  + WjCose 
o 

Normal  Tape  Equilibrium 


F^|—  = -(WjSine  + 2ocsine) 
o 


Meridional  Film  Equilibrium 

d(r^m)  = r0Wfcose  - ^ sine  [sinj  cose  + cose  sine  (cos^-  - 1)] 


[23] 


[24] 


[25] 


Normal  Film  Equilibrium 

rde_  \ _ „ n / 1 . \N  „._ir  _i_  N 

5o 

0 o 

+ sine(cos^cos':e  + sin  e)] 


Vm'ar’  ' -ro“fs1n0  - Pr<1+em^  s1"ff  * C0SB  cose 


[26] 


It  may  be  noted  that  the  governing  equations  in  the  normal  direction,  equations 
[24]  and  [26],  are  not  influenced  by  the  shear  stress.  However,  at  this  point  no 
distinction  has  been  made  between  the  angle  e when  referred  to  the  tape  or  the 
film.  A discussion  of  this  critical  feature  will  be  deferred  to  the  next  section 
of  this  report. 

Several  additional  equations  are  required  to  form  a complete  set.  If  the 
circumferential  strain  is  defined  as  the  change  in  gore  width  per  unit  original 
gore  width,  it  may  be  expressed  as: 

[27] 


1+P  - r slmr/N  (g+B) 
c r0sin(q+eT  */N 


This  equation  is  obtained  by  considering  Figure  3.  The  deformed  gore  width 
as  shown  here  is: 

* = 2V 

where  R£  is  the  circumferential  radius  of  curvature  and  y =a+g.  The  undeformed 
gore  width  is  obtained  from  the  manufacturedshape,  i.e., 

lo-Tr° 

where  rQ  is  the  design  value  of  balloon  radius  at  any  point. 
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The  angle  a is  obtained  by  noting  that  the  distance  between  load  tapes  in  the 
horizontal  plane  is  identical  to  the  distance  between  load  tapes  in  a plane 


containing  the  circumferential  radius  of  curvature,  Rc-  Then: 

-- v-  si  not  = rsin£ 
cose  N 

or  a = sin-1[sinjSjCOse]  [2f 

By  the  same  logic: 

Rcsiny  = rsin^  [ 2C. 

Equation  [27]  may  then  be  obtained  from  the  definition  of  circumferential  strain 


The  material  is  assumed  to  be  linearly  elastic  and  orthotropic  so  that  Equations 
[12]  and  [13]  are  still  applicable.  Finally,  the  geometric  relationships  are 
employed  such  that: 

35-=  0+eJsine  [3( 


= (Unloose 

It  is  now  a matter  of  arranging  these  equations  intoa  suitable  form  which 
will  be  tractable  for  solution. 


VI.  ANALYSIS  PROGRAM  DEVELOPMENT 


In  order  to  solve  the  set  of  equations  just  developed  it  is  necessary  to 
make  some  additional  assumptions  regarding  the  stress  field  and  the  tape  angle. 
Since  the  tape  force,  F,  appears  in  both  tape  equilibrium  equations  it  is  necessary 
to  employ  some  sort  of  relationship  between  this  force  and  the  film  forces. 
Therefore,  it  is  assumed  that  the  tape  force  is  proportional  to  the  meridional 
strain  in  the  film.  Therefore, 


F = KTe  = O o + KTD  o 
Tm  Tmm  Tmcc 


[32] 


Although  this  assumption  is  reasonable  at  the  interface  between  the  tape  and 
film,  the  assumption  of  a uniformly  distributed  stress  across  the  gore  will  result 
in  a uniform  strain  which  will  dictate  the  tape  force.  The  alternative  to  this 
approximation  is  to  allow  the  meridional  strain  to  vary  across  the  gore.  This 
would  result  in  a set  of  partial  differential  equations  which  would  increase  the 
computational  time  considerably.  Another  approach  would  be  to  assume  a strain 
distribution  across  the  gore  and  use  the  integral  of  the  resulting  stress 

distribution.  An  approach  similar  to  this  was  used  by  Alexander  (6). 

The  equations  of  equilibrium  may  be  arranged  in  a more  familiar  form  if  the 

total  meridional  load  is  defined  as  the  sum  of  the  tape  forces  and  the  film 
forces.  Specifically, 

NF 


r = r a + 5 — 

o m Ztt 


[33] 


The  two  meridional  equilibrium  equations  [23]  and  [25]  may  be  added  together 
to  finally  yield, 

dT  _ ,NWt  . _ „ . N 


dl~  = + r0Wf)cose  + ^csine[sin(j^)cose  + sinBcose(cos*  - 1)] 


[34] 


'o  - - ’ - - " N 

It  may  be  noted  that  the  shear  stress  which  appeared  in  both  of  the  contributing 
equations  is  self  equilibrating  and  does  not  appear  in  equation  [34].  In  addition 
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if  the  circumferential  stress  is  assumed  to  be  zero,  this  equation  will  reduce 
to  the  design  equation  for  a natural  shape  balloon  as  reported  by  Smalley. 


I 

t In  a similar  manner  the  equilibrium  equations  in  the  normal  direction, 

equations  [24]  and  [26]  may  be  combined.  However,  in  doing  so  it  is  necessary 
to  assume  that  the  tape  and  film  are  deployed  at  the  same  angle.  When  added 

. 

N 

together  the  following  equation  is  obtained: 

Iff"  = ^ + roWf)sin0-pr(1+Em)^inR  + ^c[sin=cos6cose+ 
o 

sin&(cosj|jCOs2e  - cos2e)]  [35] 

This  equation  will  also  reduce  to  the  familiar  design  equation  if  there  is  no 
circumferential  stress.  It  must  be  observed  that  in  order  to  obtain  this 
equation  no  distinction  may  be  made  between  the  tape  and  film  angles.  This 
approximation  is  quite  reasonable  over  the  majority  of  the  gore  length;  however, 
near  the  ends  e must  be  thought  of  as  the  average  of  two  different  angles. 


A third  useful  relation  may  be  obtained  if  it  is  assumed  that  the  derivatives 
in  equations  [24]  and  [26]  are  equal.  If  the  difference  between  these  two 
derivaties  is  set  equal  to  zero  then  the  following  algebraic  equation  will 
result: 

. ir  N r • " ^ / * 2„ 


rQWf  sine  + pr(l+em)^sin^  - ^°c[sin^  cose  cose  + sine(cos^  cos^e  + sin^e)] 


-57  (“t51"6  + 2»cs1"6'  - 0 


[36] 


This  equation  places  a very  severe  limitation  on  the  shape  of  the  deployed  balloon 
since  it  effectively  requires  the  tape  angle  to  be  identical  to  the  film  angle. 
Although  this  assumption  appears  to  be  reasonably  valid  over  most  of  the  gore 
length,  it  is  obviously  in  error  near  the  end  points.  However,  since  this  is 
an  algebraic  equation  rather  than  a differential  equation,  it  may  be  used  to 
determine  the  lobing  angle  for  any  set  of  shape  variable  without  integration. 
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This  equation  has  no  counterpart  in  the  design  procedure  since  no  distinction 
is  made  between  the  tape  and  film  angles. 


Since  equations  [34],  [35],  and  [36]  do  not  involve  the  shear  stress, 
this  variable  need  not  be  evaluated  to  determine  the  shape.  However,  once  the 
shape  is  found  equation  [23]  may  be  used  to  evaluate  the  shear  stress  and 
ultimately  the  principal  stresses  and  tension  field  patterns. 

In  order  to  obtain  solutions  at  a variety  of  altitudes,  it  was  decided  to 
include  a volume  calculation.  This  is  necessary  to  provide  an  overall  boundary 
condition  on  the  problem.  Since  the  total  volume  is  needed  at  the  bottom  of 
the  balloon  to  establish  equilibrium  with  the  payload  it  was  decided  to  integrate 
the  governing  differential  equation  from  top  to  bottom  with  the  boundary  condition 
of  zero  volume  at  the  apex  of  the  balloon.  All  other  differential  equations  will 
be  integrated  from  bottom  to  top  so  that  for  a positive  change  in  gore  position, 
the  change  in  volume  will  be  negative.  The  governing  differential  equation  may 
be  simply  stated  as: 

3^  = (1+em)H?  = • 0+Si)Acose  C373 

In  this  case,  A is  the  area  in  a horizontal  plane  which  is  obtained  by  multiplying 
the  area  enclosed  by  a single  gore  by  the  number  of  gores.  The  area  in  question 
is  composed  of  two  parts;  a)  the  triangle  formed  by  straight  lines  in  the 
horizontal  plane  connecting  the  two  adjacent  load  tapes  and  the  centerline;  and 
b)  the  projection  of  the  lobe  onto  the  horizontal  plane.  This  second  area  will 
be  positive  or  negative  depending  oh  whether  or  not  the  pressure  differential  is 
positive  or  negative.  It  may  be  shown  by  considering  Figure  3 that  the  area 


of  a single  gore  in  the  horizontal  plane  is  given  by: 
„2  £ n sin2(£)[2(a+e)-sin(2a+2e)] 


A = -^sin  (jj)  + cose 


sin  (a+e) 


The  choice  of  the  proper  sign  to  be  used  in  this  equation  is  made  by  dividing 
the  pressure  by  the  absolute  value  of  the  pressure. 
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Equations  [30]  through  [37]  must  be  solved  simultaneously  in  order  to  insure 
a compatible  shape.  An  attempt  was  made  to  solve  this  set  of  equations  utilizing 
the  same  modified  Runge-Kutta  numerical  technique  which  is  routinely  used  for 
the  design  of  balloon  shapes.  However,  this  technique  is  best  suited  for  initial 
value  problems  rather  than  boundary  value  problems.  This  so  called  "shooting" 
method  is  unable  to  impose  a final  boundary  condition  and  must  rely  on  the 
proper  adjustment  of  initial  conditions  until  the  desired  solution  is  obtained 
at  the  final  boundary.  In  this  particular  case,  the  solution  of  five  coupled 
differential  equations  is  stable  for  only  certain  values  of  circumferential 
stress.  As  sufficiently  large  stresses  are  developed  the  integration  scheme 
becomes  unstable,  diverges  and  a solution  is  unobtainable.  It  is  for  this 
reason  that  limits  are  placed  on  designs  with  circumferential  stress.  It  must 
be  emphasized  that  this  is  a difficulty  with  the  numerical  technique  and  not 
with  the  formulation  of  the  problem. 

In  an  effort  to  obtain  a stable  solution  with  implicit  control  of  the  boundary 
conditions,  the  governing  differential  equations  were  rewritten  in  finite  difference 
form  to  form  a set  of  nonlinear  algebraic  equations.  The  balloon  gore  is 
divided  into  a large  number  of  points  and  the  six  equations  written  at  each  point. 
This  results  in  a set  of  equations  which  can  be  efficiently  handled  due  to  the 
banded  nature  of  the  problem.  A Newton-Raphson  technique  was  then  attempted 
but  did  not  yield  a convergent  solution.  It  is  difficult  to  ascertain  the  reason 
for  this  lack  of  convergence  but  it  is  suspected  that  the  derivatives  of  the  non- 
linear function  defined  by  equation  [31]  near  the  top  of  the  balloon  approach 
zero  which  cause  the  coefficient  matrix  to  become  singular. 

A direct  iteration  technique  was  employed  to  eliminate  the  instability 
associated  with  the  inversion  of  a singular  matrix.  In  this  technique  an 
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equation  is  written  for  each  variable  in  terms  of  known  quantities.  For  example, 
equation  [30]  may  be  used  to  evaluate  the  radius  at  point  i,  r^ , for  the  next 
iteration  step  according  to  the  equation: 

ri  ■ vi  • 4s„o%>is,nei 

Similar  equations  are  written  for  each  variable  at  every  point.  The  right  hand 
side  of  each  equation  is  evaluated  from  current  values  and  the  left  hand  side 
is  used  in  the  next  iteration  step  on  the  right  hand  side  until  the  changes  are 
suitably  small.  As  in  the  case  of  the  Newton -Raphson  technique,  the  direct 
iteration  technique  was  found  to  be  divergent  when  applied  directly.  However, 
by  using  an  under  relaxation  technique  where  only  a portion  of  the  changes  called 
for  in  each  variable  were  used,  the  solution  became  stable  and  convergent.  When 
the  design  shape  is  used  for  the  initial  shape  estimate,  convergence  is  obtained 
in  13  iterations.  An  intermediate  step  in  the  iteration  process  is  the  computation 
of  meridional  and  circumferential  film  stresses  at  each  point  along  the  gore.  This 
is  accomplished  by  solving  equations  [32]  and  [33]  simultaneously  with  the  membrane 
equation  applied  at  the  center  of  each  gore,  i.e., 

^ + pr-  * p + Wfsine  [39] 

The  angle  e is  positive  when  measured  clockwise  from  the  vertical  parallel 
to  the  balloon  centerline.  As  a result,  since  s is  measured  positive  from  the 
nadir,  the  meridional  radius  of  curvature  will  be  positive  if  the  change  is 
angle  is  negative.  Therefore: 

1 _ de  1 
" dso 

In  equation  [39]  the  weight  per  unit  original  area  is  assumed  to  be 
negligibly  different  from  the  weight  per  unit  actual  area. 


VII.  ANALYSIS  PROGRAM  RESULTS 


The  analysis  program  was  applied  to  the  balloon  described  in  Section  IV 
of  this  report.  Appendix  B contains  a detailed  description  of  the  initial  shape 
and  weight  distributions  as  well  as  the  deformed  parameters  given  by  the  usual 
design  process.  This  particular  balloon  was  selected  for  analysis  because  it 
is  very  similar  to  the  balloon  analyzed  by  Alexander  (6)  in  attempting  to 
determine  the  circumferential  stresses  due  to  lobing.  In  addition,  this  bal- 
loon is  typical  of  the  "heavy"  load  balloons  which  have  experienced  an  inordinate 
number  of  failures  in  recent  months.  The  problem  was  first  solved  for  the  design 
altitude  so  that  the  results  could  be  compared  to  those  of  Alexander.  It  was 
then  assumed  that  the  balloon  was  in  equilibrium  at  a variety  of  altitudes 
both  above  and  below  the  design  altitude.  The  altitude  regime  below  the  design 
altitude  was  selected  to  correspond  to  those  altitudes  that  would  have  been 
recommended  by  the  manufacturer  on  the  basis  of  cap  thickness  and  load  tape 
strength.  A convergent  solution  was  obtained  for  each  of  the  altitudes  of 
interest. 

At  the  design  altitude,  the  solution  converged  to  a payload  ratio  of  1.001. 
This  indicates  that  the  effect  of  lobing  is  to  increase  the  volume  slightly 
over  the  design  volume.  The  lobing  angle,  8,  varies  continuously  along  the 
gore  but  reaches  a local  minimum  at  a point  between  the  maximum  radius  and 
the  edge  of  the  cap.  The  maximum  amount  of  lobing  occurs  at  the  top  of  the 
balloon.  The  departures  from  the  design  values  of  radius,  height  and  volume 
are  minimal  and  as  a result,  all  computations  must  be  performed  in  double  preci- 
sion. The  most  significant  departure  is  the  angle,  0,  and  the  total  load  near 
the  bottom  of  the  balloon 

The  meridional  and  circumferential  stress  distributions  are  shown  in 
Figure  6.  The  results  reported  by  Alexander  are  also  shown  in  Figure  6 for  com- 
parison purposes.  It  should  be  noted  that  the  current  results  have  been 


dimensional i zed  with  respect  to  the  shell  thickness.  This  assumes  that  the 
cap  is  not  a load  carrying  member  at  the  design  altitude.  It  is  possible  that 
the  large  discontinuity  in  Alexander's  results  at  the  edge  of  the  cap  is  due 
to  not  only  the  weight  discontinuity  but  his  assumption  that  the  cap  is  a 
load  carrying  member.  However,  Alexander's  results  indicate  that  the  circum- 
ferential stress  increases  as  one  approaches  the  apex  of  the  balloon.  It  is 
felt  that  the  present  results  are  more  reasonable  since  the  circumferential 
radius  of  curvature  is  proportional  to  the  radius,  r,  and  should  approach 
zero  as  the  apex  is  approached. 

The  results  of  the  present  formulation  are  presented  in  Figure  7 for 
a variety  of  altitudes  of  interest.  The  meridional  stresses  are  shown  only 
for  b/bp  = .8  and  b/bp  = 1.4.  At  any  intermediate  altitude,  the  meridional  stress 
will  have  a value  between  these  two  curves.  This  narrow  band  is  relatively 
unaffected  by  changes  in  altitude  of  interest.  However,  the  circumferential 
stress  distributions  at  four  different  altitudes  are  also  shown  in  Figure  7. 

It  should  be  noted  that  at  altitudes  below  the  design  altitude  the  circum- 
ferential stresses  increase  significantly  above  those  predicted  at  the  design 
altitude.  Above  the  design  altitude,  in  a regime  of  decreased  pressure  differ- 
ential, the  circumferential  stress  is  significantly  reduced. 

Although  the  stresses  predicted  by  the  present  analysis  are  in  themselves 
not  considered  severe  enough  to  cause  failure  of  this  balloon,  the  region  in 
the  shell  below  the  edge  of  the  cap  is  a region  of  large  biaxial  stresses. 

Any  flaw,  damage  or  structural  discontinuity  in  this  region  could  cause  a 
sufficient  amplification  of  these  stresses  to  result  in  uncontrolled  deformation 
or  even  failure.  This  region  will  be  most  severe  at  different  times  depending 
on  the  time  of  launch.  If  a morning  launch  was  attempted,  the  highest  magnitude 
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Figure  7 - Stress  Distribution  in  Off  Design  Configurations 


of  stress  will  occur  as  the  balloon  with  all  ballast  still  on  board  goes  to 
its  equilibrium  altitude.  As  ballast  is  dropped,  the  balloon  will  rise  and 
the  stresses  will  decrease. 


Should  this  balloon  be  launched  in  the  evening,  the  system  may  assume  an 
equilibrium  altitude  with  a subpressure  region  due  to  the  low  gas  temperature. 
However,  at  sunrise  the  gas  will  expand  to  fill  the  available  volume  and  the 

pj 

balloon  will  rise  to  assume  the  maximum  stress  condition. 

l 

1 


VIII.  CONCLUSIONS 


The  problem  of  design  and  analysis  of  a single  cell  balloon  has  been 
reformulated  in  a manner  which  will  permit  a more  realistic  determination  of 
the  state  of  stress  in  the  thin  film.  The  design  procedure  has  been  modified 
to  yield  a first  order  estimate  of  the  meridional  stress.  However,  the 

I 

stresses  computed  in  this  manner  are  not  conservative  and  second  order  effects 
must  be  considered  . The  effects  of  lobing  have  been  included  in  an  analysis 
procedure  which  produces  a realistic  distribution  of  both  meridional  and 
circumferential  stress.  The  model  formulated  in  this  report  suggests  a 
plausible  explanation  for  the  transfer  of  loads  between  the  load  tapes  and 
the  film. 

A computational  technique  has  been  developed  which  is  capable  of  solving 
the  equations  formulated  in  th's  report.  The  equations  are  highly  nonlinear 
and  no  generalizations  may  be  made  nor  can  the  uniqueness  of  the  solution  be 
guaranteed.  However,  results  have  been  presented  for  a typical  heavy  load 
balloon  of  interest  to  many  organizations.  It  has  been  shown  that  the  cir- 
cumferential stress  is  significant,  may  exceed  the  meridional  stress  under 
certain  conditions,  and  achieves  its  maximum  value  at  the  edge  of  the  cap. 

It  is  hoped  that  the  technique  presented  in  this  report  will  find  accept- 
ance by  those  interested  in  the  successful  flight  of  high  altitude  balloons. 

Many  observed  features  such  as  slack  load  tapes,  stress  bands  and  lobing  may 
now  become  predictable  events. 
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LEVEL  *1 


MAIN 


DATL  = 7G  1 0 1 


1 


C* *********************  **************************  A*************  ***** ♦ ** 

C - b A L D L - 

C BALLOON  DESIGN  PROGRAM 

C*  •• *******  ********  ************  * * *******  *********  ******  *******  *********  * 

L IMIS  PROGRAM  CALCULATES  THt  SHAPE  OE  A FREL  UALUUUN  GIVEN  A Si  1 Of 

C INPUT  CONDITIONS.  IT  IS  EUR  EULLV  TAlLURLO  BALLOONS  ONLY.  W 1 T H 

C WITHOUT  AN  ENDCAP.  C ONVERGLNCt  UE  TEC  SOLUTION  IS  ObTAlNtU  bY 

C MATCHING  TMt  IUP  L*jAL>  REACH  UN  TU  T ML  TOP  LOAD.  IMIS  PRUGRAM 

C INCLUDES  TMt  EFFECTS  UE*  1 Mt  LOAD  I APES  IN  TME  CALCULATIONS. 

C* *********  ************  ************************************************* 

C THIS  PROGRAM  IS  OASt O UN  1 Ml  PROGRAMS  "WALLOON"  AND  "ELATLALC". 

C AHCRL-CS-92.  "DLILWM I NATION  OF  TME  SHAPE  UF  A FREE  BALLOON".  HV 

C J.  M.  SMALLEY. 

q************** ***********************************************  ********** 

C DESCRIPTION  UE  DATA  OtCK: 

C * EACH  SE  1 UE  GAIA  CuNSISlS  UE  THKlt  CARDS  LISTING  T ML  VALUES 

C IN  T Mt  READ  STATEMENTS  SOO  AND  SOI 

C ALL  DATA  ARL  REAL  IN  10  wluL  FIELDS  OF  TMt  APPROPRIATE  TVPt. 

C 1 "FIRST"  LARD  INPUTS  AHt  AS  FOLLOWS ! 

C PAYLOAD  IN  PuONLS 

C ALT  1 1 UD L OPTIONS  1 - ALT  IN  FT,  2 = ALT  IN  MB 

L ALTITUDL  IN  FEET  UR  MILLIBARS  AS  DESIRED 

C E 1LM  lYPL:  1 = PULYLTMYLENL  • 2 - MYLAR 

C FILM  THICKNESS  IN  INCHES 

C LOAD  TAPE  TYPES  1 - POLYESTLR,  i.  = KEVLAR 

C I APE  LUAO  RATING  IN  POUNDS 

C NUMUE  R LF  LUAO  TAPES 

C 2 "SECUND"  CARD  INPUTS  ARE  AS  FOLLOWS: 

C TOP  LUAD  IN  POUNDS,  ( ■♦  ) UP.  (-)  DOWN 

C STRESS  CONSTANT  TaUO  l USUALLY  0.001 

C STRESS  CUNSTANT  TAU1  (USUALLY  0.00) 

L SUPtKI’Rt  SSoRt  10.00  FOR  NATURAL  SMAPL  ) 

C PRINT  1NLRLMENI  N (0  FOR  S1AMIARD  OUTPUT) 

L NUN— O IMLNSIUNAL  GURL  INCRLMLNT  DSO 

C NON-DIMENSIONAL  GURE  LENGTH  TCI  CAP  STARTING  LOCATION  CSTAUT 

C FILM  THICKNESS  INCLUDING  LAI'  IN  1NCHLS 

C J "THIRD"  LAkD  INPUTS  AHt  AS  FOLLOWS: 

C OUTPUT  CONTROL  (K!  V^-2  FuR  PUNCHED  DLCK  OF  SHAPE  L Wt  loHl 

C KtV2=l  FOR  DISK  FILC  OF  SHAPE  & WLlGHT) 

C OUTPUT  CONTROL.  (MPT  IS  NUMOtR  UE  POINTS  IN  LOAD-  AL I I I Uofc  » 

C IDENTIFY  LIFTING  gAS  (FUR  MCL1UM,  1GAS-1I 

C MINIMUM  UtCUMMt NDL D PAYLOAD  IN  POUNDS 

C MAXIMUM  KECOMMENOI  D PAYLOAD  IN  POUNDS 

C * IE  NO  ENDCAP,  MAKE  SURE  THAT  CS TART  IS  GREATER  THAN  ANY 

c anticipated  gore:  llngtm. 

C * TMt  PROGRAM  REOUIRt-.S  A LAST  DATA  SET  WITH  P = 0.0  TO 

C TERMINATE 

C**A  ***  *****  ***************** **♦ * **♦*•♦*  *****  **********  ***************** 

C DC  SCR  I P 1 1 ON  OF  OUTPUT: 

C FOP  EACH  CASE  TME.  INPUT  OAT  A WILL  FIRST  HE  PRINTED.  FOLLDWLO 

C BY  A RECORD  OE  TME  IIEKAT1UNS  R LOU I RED  FOR  CONVERGENCE  OF  TMl 

C SOLUTION.  THIS  RECORD  CONSISTS  OF  A DISPLAY  OF  INITIAL  AND 

C FINAL  ANGLES  OE  THE  GOWt  WRT  THE  VERTICAL  AXIS  OF  TME  BALLOON. 

L DURING  TME  FINAL  ITERATION,  PERTINENT  VALUES  ALONG  Till.  Gt-Kfc  A~ 

C PRINTED  OUT  CUNT  I NOUUSL  Y . FOLLOWED  BY  A LISTING  DF  FINAL  VALU: 

C***  ******************************************************************** 

C APPRUPRIATc  ERROR  Mt SSAGtS  ARC  SUPPLIED  IF  THE  SOLUTION  DUL S NOT 

C LUNVtRCL 
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v r 


r\Cr i n no 


tVtL  21 


MAIN 


DA  It 


76226 


2 l /?'»/42 


******  **************** 


MAIN  PROGRAM 


RtAL  K .KT .KT APE  .L AMEDA . I AM 
iNTtGER  C COE F »C  CDE T 

DIMtNS ICN  GP(100).G*(100).CWT( lOO).GT(lOC) 

100  WRITE(6.6C0> 

600  FORMAT! • 1 • ) 

RE AD  I 6 .5  0 0 > P . K£ Y , C C NS T , CCCE F , FT h l CK , CO DL T , T LP . N T 

600  FORMAT(F1C.0.110.F10.0.I10.F10.4,I10.F1U.0.110) 

READ! 6 • 50 1 > TL. T ADO. T AU 1, ALP  HA, N.DSO.C START, LAP 

601  FCFMAT(F10.2.F10.2.F10.2.F10.2.110,F10.2.F10.2.F10.4) 


IF  NC  MORE  DATA,  EXIT. 
1 F ( P ) 2CC.3CC.2C0 


200  CCNT1MJE 

R t AO  ( 6 . 5 0 2 ) KL  Y 2 . VFT . IGAS.FMIN.PMAX 
602  FCPMAT I 3 1 10 . 2F 1 0 .0  ) 

CALL  CE  SIGN  (F.KEY.CCNST.CCCEF  .FTHICK.  CuoET  .TLR.  NT  , IU  . UUO,  T AJ  1 , 
* ALP  FA ,N.DS0.CSTART,CAP,GP.G*,WbAL,KtY2.S) 

CCMFLTE  LCAC  - ALTITUDE  C J BV E 

wRITt (6 .600 > 

WR  I Tt  ( 6 ,6C1  ) 

601  FCRMAT150X.  'LOAD  - ALTITUDE  D A TA  • , / 4 2 X , * GRUS S A l POOR N*  */37X»* 

I'aEIGhT  - (KG)' .10X.' ALTITUDE  - (KM)') 

GPMIN=l.-M»dAL 
GPMAX=GFMIN-MPMAX-PM  IN  )/P 
OG=(GPMAX-GPMIN ) /MPT 
CP  1=CPM  IN 
260  CCNTINUE 

CALL  GLNoTHCGPl  .SLMDA.l) 

L’1=P/ (S/SLMDA  )*  *3/.  0 62  43 
CALL  DCYNCY ( El  .P A . I CAS .T CFIGH ) 

CALL  ATMGS2(PA,H,TA.RHCA.e.lAM,GXPAN.T0HlGP> 

GF  CS  = CP  1 *F*  .464 
H=H/1  CCC  . 

WRITE (6.602  ) GR  OS.H 

602  FCRMAT(42X  .E15 .7  .8X.E15.7  ) 

GP 1 =GP1 FOu 

1F(«P1  .GT  .GPMAX  )G0  TO  100 
GC  TC  250 
300  CONTINUE 
STLP 
END 
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NHWMf 


or'C  o r nn 


EVEL  21  DESIGN  OATt  = 78226  21/29/42 

SLBRCLTINE  CESIGN  (P.KEY.CCNST  .C COEF  , FT H I CK , CCDcT , TL R. NT , TL * TAUO, 
STAUl.ALPHA.N.OSC.CSTART.CAP.GP.GW.WBAL.KEYE.S) 

**********  *****1 if*  **************************************  **************** 

_******************  SOORCUTINE  DESIGN  **************  * *********  * 

REAL  K.KT  .KTAFE  , NU 
WEAL  LAMBDA 
INTEGER  CCDEF.COOtT 

D1 MENSI CN  Y(6).C(6).G(6).E(4),f(4).D(4 ),WCRD1 ( 3 > , WORD2 ( 1) 

C 1MENS  ION  GP<  1 00  ) .G«(  1 CO) .GwT(  100)  .GT(IOG) 

DIMENSION  SSTCRE  (3 JO ) . CSTCRE < 30  0 >. WSTORE! 3 00  ). T STORE!  300 ) 

WEAL  Fit  3)/'PCLY'  « ' E TH  Y • *•  LENE'  / 

REAL  F2( 3 )/  'MYL A • , 'W  ','  •/ 

REAL  T1 <3)/'  FCLY*  .' ESTE'  . *R  •/ 

REAL  T2< 3) / 'KEVL • . • AR  ','  •/ 

REAL  Kl/'  FT'/ 

REAL  K2  / ' MB'  / 

d.  C.  AND  L>  ARE  CONSTANTS  USED  IN  GILL'S  MODIFIED  R-K  METHOD  FOR 
SOLVING  THE  GOVERNING  DIFFERENTIAL  ECUAT1CNS. 

DATA  6/2  .Of  1 .0.  1 .0. 2.0/.  C / 0 . S • 0 . 292 89 . 1 . 7 C7i  . C. 5 / . D/0.6. 

SO. 2 92 89. 1 .7071  ,0 .16667/ 


CALL  MATL(P»KEY,C  ON  ST . CODEF • FTH1 CK.COOET  ,TLW  .NT  .CAP.SIGMA.CSIGWA  , 
STSIGMA  ,KT  ,£M ,EC .EMC.LAMUDA.TS  I GM E ) 

NU=£MC/EM 

QMEGA1=S  IGMA/l  .8453 
OMEGA2=CMEGAl*(C SIGMA/ SIGMA) 

OMEGA=OMEGAl FTS IGMA/l ,e453 
TCMLOa  = TS  I G M A/6  .2832 

EMPIRICAL  CURVE  FITS  FROM  SMALLEY'S  REPORT.  TCI  GE  T oUOD  FIRST  ESTIMAT 
UF  THETAO.  GCCD  FOR  SIGMA  BETWEEN  0.0  ANC  0.8  . 

ThETmC=I.347 

IF(LMEGA.LT  .0.43354)  T FE  T A 0*  0.72  *OME  GA+1.1015 
1 F ( C MEGA  . L T . C . 3 5225  ) TFETA0  = 1 . I 5 * CMC  G A«-o  . 950  7 
I F ( uM  EGA  .LT.0.27CS6)  THE  TA 0= 1 . 46*0M£ G AF 0 . et8  7 
(JMLGA=CWEGA1 

I0U1T  = 0 
LSTR0N=1 

lF<«_nuLF.EQ.2)  GC  TC  52 
DC  51  1-1.3 

5 1 WCRJI  < 1 ) =F 1 ( X ) 

GO  TO  54 

52  DC  53  1=1.3 

£3  WORD 1 ( 1 ) =F2<  I ) 

54  I F ( CODET  .EG  .2  ) GO  TO  5t 

DC  5b  1 = 1 .3 

55  WORD  2 < I ) — T 1 < I ) 

GC  TC  58 

56  DC  57  1=1.1 

57  WORD  2 ( I ) = T 2 ( I ) 

58  IF  (KEY.  EC. 2)  GO  TO  5‘i 
W0RD3=K1 

GC  TL  6 C 

59  WCRDJ =K2 
tO  CONTINUE 
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eVEL  21  CtSIGN  L>A  Tt  = 7b226  21/2P/42 

WR  I TF  < 6 .6  Cl  ) ( wCROl  (1)  .1  = 1.31  .FTHICK.SIbKA.TA.JO.  (WDR02  ( 1 ),  1 = 1.3). 

*EM.  T AUl  . TLH  .eC.  ALPHA 

oOl  FCPWAT(5X..»7  ( • * • )/5X , • , 35X  • . 1 5X  . 'NDN-D  I MENS  IL’NAL  Q'JAN  T I TIES*. 

lHX.'MISC.  PAKAMETEKS‘/5X.‘*  FILM  TYPE'.HX.':  • , A4  , A4  . A4  , ' *‘/SX, 

2**  F ILK  THICKNESS  J • .F  7.4  , • IN  ' .4X,  • * ' , 1 7X.  'SIGMA  FILM  =‘.FR. 4 
S.lk'X.'T  AUC  = • »Ft>  . 3 /5X  * TAPI.  TY  PE  • ,bX  • • S 1 , A4.  AA.  A A « 1 * • . 1 7X  , 

A'feM'.SX.  *=‘»Fd.3»13X»  • TAU  I = • .F  6.  3/EX.  • A TAPE  LOAD  RATING  : •, 

5F  6*0,'  LBS  *•  . I 7X . • EC  ' .RX • • FB. J, 1 jX.  • ALPHA  = '.F6.3) 

WRITE (6.602)  NT.fcVC.USC.F.TSIGMA.CSTART.TL .KT ,N, CONST . WCR03 
oO  <:  FORMAT  ( * 1 « AX  • ' * NUMEfcR  UF  TAPES  J '.IS.  OX  . • * • . 1 7 X • *F.  KC  • • PX  , 

1 ' =•  .F  6.  3 .137  .' DEO  = • ,F6 . 3/SX  . • * PAY  LOAD  '.  10X  ,'  : '.F6.0,'  LHS 

2 * • . 1 ?X.  • S 1GMA  TAPE  = • .F  S . A , l 2 X . »C  S T AH  T = • ,F  S . 2 /S  X . • * ( 1 NCI.  TOPLC.A 

3J  CF ) J • ,Fu .0 . * LBS  * • . 1 7X.  'KT '.9X.  •= • . F6 .1 . 15X.  »N  = '.I2/jX 

A,  'A  DESIGN  ALTITUDE  i '.F7.0.A4.'  * • /S  A . ' * * , 35X  . • * • / 5X  . 3 7 ( • <•  • > 

5.///SX,  ‘ITERATION  HECORC  (INITIAL  ANO  FINAL  ANGLES)*//) 


IF  THE  SOLD  T I UN  IE  NOT  CONVERGING,  EXIT 

21  I CO  1 T = I Co  1 1 ♦ 1 

IF( IG01 T.LT.20)  GC  TC  22 
wR  1TL ( 6, e03  ) 

603  FCFMATI'-','  ***  SJLJTICN  IS  NOT  CONVERGING  - MAXIMUM  NUM3FR  OF  IT 
SERAT1GNS  HAS  BEEN  EXCEEDED  ***  ») 

GU  TO  5 

LSTRUN  = 0 IS  THE  KEY  FOR  LAST  ITERATION.  IT  IS  AN  OUTPUT  CCNTRCL. 

22  IF (LSTRUN. NE  .0 ) GO  TO  2J 


*R ITE ( 6. 604  ) 

fcOA  FORMAT ( *0* »9SX  » 2 7F  M ANU  FACTURER  *S  T AUL  E L A YUUT // I 0 X . • S • • 1 3X , • P • . 

1 1 2X,  • l'  ,13X.'T‘  .llX.'TAUK*  .OX  . • T AUC*  .8X. • GCRE  POSITION  ( FT  ) • .OX, 

2 ' R AL  F CORE  WIDTH  (IN)*/) 

23  DC  24  1=1  ,o 

24  Q(  I ) = 0. 

231  JU  241  1=1,300 

SET  OR  E(  I ) = 0 . C 
W S TOR  E ( I ) = 0 . 

T S T CR  t ( I ) = 0 . 

241  G STOKE ( I ) =0. C 
J COUN  T = 1 
OS  = cso 
SSTOP  = C ST  ART 
OMfcGA  = CMtC-Al 
TCMEGA=TSIGKA/6 .2832 
T w=  0 . 

IT AG=0 
RC  = C. 

GW  e—  0 • 

SO  = 0 . 

SOU  I M = 0 » 

R = 0 . 

S = 0. 

7=0. 

WSU  V=0 • 

A S L v=  C . 

T AUM=C. 

T ALC=  TAUO 
THE  TA  =THE  TA  C 

RT0=1  ./6  .28J2/C0b<  THETA) 
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eVEL  21  DESIGN  DATE  = 78226  21/24/42 

V ( 1 >=THETA 
Y (2) =R 

Y( 2 ) =2* ALPHA 

Y (4)  = l(P*TL)/P)/t«2 832 /CO  S( T HETA  ) 

Y ( £»  = C. C 

Y(6)=C.C 

SSTCREt JCCUNT  ) = SOCIM 
GST  Of-'  E ( JCCGNT)  = G*2 
2 T =6  .2832*Y(  4 > 

IF(LSThUM34  « $2  .34 

32  wW(TE(fc.6C5)  S . R . 2 • T . T A UM  , TA  UC  . S 03 1 M ,G  W2 
605  F U A M A T (4(4X.F10  . b).4X,F3.3, 5 A * F B . 3 . 1 2 A » F 7.  1.  1 7X  »F  7.2) 
IF(KFY2.EQ.l)«RlTE(3.653>Y(2>.Y(l),Y(4).Y(l).Y<6) 

1F(NEY2.Eu.2  ) WRITE!  T.6  50)Y<2).YJ2)  ,YI4)  ,Y(1)  ,Y(6) 
o50  F CRMAT(5£15  .7 > 


1CCUM  AND  N ARE  OUTPUT  CONTROLS  USED  DURING  THE  LAST  ITERATION. 

N IS  THE  NUMBER  CF  INCREMENTS  THAT  « ILL  BE  SKlPPEO  IN  THE  PRINTOUT. 

34  ICCLNT=N 


IT  AG  SENDS  THt  PROGRAM  INTO  ITERATIONS  AND  RELEASES  IT  WHEN  THE  TOP 
CF  THE  BALLCCN  HAS  EEEN  REACHEC 
I F ( 1 T AG- 1 ) 4 . 3 . 4 


3 I F I LS TRUN.EC . C 1 GC  TC  33 
OEGO  = THETAt*S7.29Sc 
THETAL  = THET A*57 .2956 
DS  = DSC 

»R  1TF (6.61 1 ) OEGO.THLTAD 
611  FCRMATC  • ,5X  .F  10  .5 ,5X  .FI  0 .5  » 


CONVERGENCE  TESTS  TOP  LOAC  REACTION  WITHIN  .5  LH  OF  APPLIED  TCP  LOAD. 
F = -e.2P32*Y(41*C0S< THETA)*P 
IF(  ABSI  F*TL  ) .LE  .0.50  LS  TRLN  = 0 

CCMPGTE  CORRECTIVE  TERM  CN  THET  AO 

DELTAITHETAO  l/DEL  TA  ( THETA!  IS  GENERALLY  APPROX. = (-F1/-2) 

D THET  A IS  THE  CESIREC  VALUE  LF  THETA  TU  MATCH  THE  TUP  LOAO  RCACTION 
TO  ThE  APPLIED  TCP  LOAD. 

CTHETA  = -(  ARCOS(TL/(6.2832*P*Y(  4)1)  ) 

DELTA  = DTHET  A - THETA 
CCRR  = -DELTA/2. 

T HE  T AO  = THETA04C0RR 

GO  TO  2 1 


COMPUTE  AND  PRINT  FINAL  QUANTITIES 
35  VOL  = J.1416*Y (o 1*(LAMEDA**3 ) 

ASUM  = ( ASUMHE. 2632*Y(  5)  )*(LAMBDA**2  ) 

WSUM  = ( W SU  M 46 ,2832*UMEGA*Y{ 5)  ) *P 

TW=TW*6.2632*P 

WEIGHT  = iSUMtTm 

F = -6.28324Y (4 ) *COS (THET A ) *P 

THE  TA 0 = 57. 2956*  THETA 

DECO  = THET A0*57 .2956 

J M AX  = JCCUNT 
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-EVtL  21  UESIGN  DA  ft.  = 7 82  2 C 'M/29/42 


PTCF=-TL/F 
*OAL  = ME  IGHT/P 

wR  IT  t (6  .o06  ) LAMBDA . VOL . DEoO  * A SUM* THE  TAD . wSUM ,f  . TW  . *E I GH  T 
t Ob  FORMAT! ,4X  . *F I NAL  QUANTITIES  (LAMBDA  =*.F7.1,*  h T ) : ' // 5 X.  • VOL  UM 
It  • .6X  • • = • «E  1 0.  3.  ' F T 3 • . 1 0 X < ' INIT  IAL  A NGLfc  • .5  X • • = • ,F  / . 3 . • OEG'/SX, 

2 • AREA  • ,8X. ' = * . tl 0 .3.  • FT2  • * 10X  , *F  INAL  ANGLL • • 7 X . • = * . F 7. 3 . • DEG* 

3/ 5 X • • FILM  HEIGHT  = ',F7.0.*  LOS • . 1 OX  . • T CF  LOAD  REACTION  ='.F7.3. 

A*  LOS  */EX.  • TAPE  HEIGHT  =',F7,Oi'  LOS* /iX.  • TOT  . WEIGHT  = *,F7.0. 

5*  LBS'/) 

IF(KEY2.E0.1)wHlTE(3.650>F,PTaP.WEAL.S.LAVuDA 
I F { KE  Y2  .EC  .2  >\*R1  Tt  ( 7,  6 5 0 ) P . P TOP  , wtlAL  . S.LAMflUA 
IFIKLY2.EC.1  ) W R I T E (3 *o51  ) NU. KT. E*.NT 
IF(KEY2.EU.2)WWITE(7.6S1)NL.KT,EW.NT 
051  FORMAT < JL15  .7,  I 15  ) 

GENERATE  TRANSFER  MATRICES  ...  ICC  FVE  NLY  SF’ACT  D . UNDtrORMCO 
GCRl.  CGCRC  INATt'S 

USNtW  = SSTUREI JMA X ) /5S. 

J=1 

SNEW=C. C 
GP ( J ) = 0 .0 
GM  J)  =0.0 
GwT!  J ) = T C ME G A 
GT ! J > =T  CMfcGA*6 .2  832/NT 

WRITE (b  .600)  J.GP!J).J.GW(J).J.GWT(J).J,GT(J) 

100  SNEW  = SNEwEOSNEw 
J = J + t 

UO  1 C 1 1=1.  JMA  X 
(TUAN  = SNLW- SSTORE  { I ) 

IF (ULAN.LT.O .0 ) GC  TC  102 

101  CONTINUE 

102  GPIJ)  = SNLw 

Gw!  J)  = uSTCRE ! I -1  > ♦ ! ! SNE w-SSTCRE!  1-1  > > + (GSTLWL!  ! >-GSTURL(  I-  ! > )/ 
t( SSTUREI  1 >-SSTORL<  I- 1 ) ) ) 

G*T ! J ) = WSTCRt  ( I - 1 ) ♦ I (SKEW -S S T OR E ( I - 1)  > * ( * STORE!  I >- w STOP t ( I - H ) / 

*(  S STORE ( I 1-SSTOHE I I -1)  >) 

GT  ( J ) = T ST  ORE  I 1-1  )■*•(<  SNE  W-SSTuR  L(  I - 1 ) )*  ! T STURE  ( I ) - TSTURl  I l - 1 ) )/ 

* I S STORE ( I ) -SSTCRE I I -1  ) ) > 
wRITfc  (6.600)  J.GP(J).J.GMJ).J.G*T(J).J,GT(J) 

600  FORMAT!*  *.5X.  • GP ! *.  13.  • ) - • . F l 0 . 3 • 1 OX . *Gw(  * • I 3 , * ) = • *F 1 0. 3 . 1 OX . 

1 * G H T ! * »I3»*  ) = • »F1  0.3  ,1  OX  . *GT  ( • ,1  3 . • ) = *.F10.3) 

652  FORMAT! 4E15 .7 ) 

IE  ! J. tO. 100  ) CC  TC  103 
GU  TC  ICO 
10  3 C C NT  1NUE 

DC  300  J = 1 . JMAX 

IF(KEY2.Eu.l)WRITE!3.652)  S ST  ORE  ! J)  .G  STORL  ! J)  . W ST  CRE  ! J ) .T  ST  UP  E ! J ) 
IF(KEy2  . EC  .2  )WRITE(7,652)SST0RL!  J ) . GSTORE  ! J)  . W STORE!  J)  . T STORE!  J) 
300  CONTINUE 
C 

GC  Tu  5 

CUMPUTAT IONAL  PORTION 
C 

C CEECK  IF  CVER  THE  ShUULOER  OF  THE  BALLOON 


4 I F ! SI  N!  THET  A ) )6  .8  ,8 


49 


C CHtCKS  FUR  IMPUSSI-JLfc  liALLCLKS. 

C IF  |TK1TAJ>PI  UR  S>10.0.  PRINT  AN  A(5°WOPfi  1 ATE  ERROR  MESSAGE  AND 

C RHCCEEU  TC  THE  NEXT  SET  CF  CATA 


-—.A 


.......  . . . . 
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EVEL  21 


JE  SIGN 


DATE 


7U226 


? 1 /?5/4  - 


17  IF<  <AGS(  THE  TAI-3.141&I  . LT.O.OI  GC  TC  12 
*R  IT  E ( 6 • 607  ) 

607  FORMAT!  ••  ***  SOL  JT  I CN  IS  NCT  CONVERGING  - IMPOSSIBLE  BAL  LOCN  : 

$ |TMETA|  > PI  ***  •* 

GC  TC  5 

12  I F ( < S-l C. C) .L T.O  .0  ) GC  TC  30 
»RITC(o.6Ce> 

6 OB  F LRMA  T<  * — * • * ♦**  SOLUTICN  IS  NOT  CONVERGING  - IMPOSSIBLE  BALLOON  : 

S NON—  O I M S > 10.0  •> 

GO  TO  S 


THIS  PORTION  ALLOTS  FOR  THE  CHANGE  IN  F1LK 

20  IF!  S-SSTCF1  14  ,13.13 
12  T=6»2632*Y!  4 ) 

131  ASUM= ASUM*6 .2832 *Y (5 > 

WSGM=  »S(JMF6.2  832*  CVEGA4Y  (5  ) 

Y ! 5 >=0  . 

I F ( LSTRUN.NE .0 ) GC  TC  132 
60*  FORMAT!  *0'.  3X.  'CAP  STARTS  AT  S = '.F6.2,  • 

WRITE  (b ,609  ) CST ART . CS IGM A 
122  G«EGA  •=  CKEGA2 

T GMEG  A=  T 5 IGMb / 6 • 2 E 3 2 
S S TCP  = 10. 0 

14  I F ( ICOUNT  >2  *2  .IS 

15  ICOUNT=  ICOUNT- 1 

II  F ( i T AG  ) 2 ,4  .2 

5 CONTINUE 
S=S*L  AMEUA 
RETURN 
t Nl) 


THICKNESS  DUE  TO  THE  CAP 


o nn  r r n nr  nr.  r nr.  n nrr  r nr>onr  or  rmi  >< . 


EVtL  2 1 


matl 


DATE 


7H22S 


2 1 / 29 / 42 


bUPROUT  1NE  MATL  <P, KEY,  CON  ST  . CO  UE  F . F T F IC  K , C.  CiU  tT  • T L R .NT. CAP,  SIGMA, 
SCSIoMA.TSIGMA.KT  • b M , t C .EMC  iL*^6CA,TS  I GMB  > 

>++ I**************************************************************** 

THIS  SUBROUTINE  CALCULATES  MATERIAL  PROPERTIES  GIVEN  I NF  l -i  MA  T I CN 

AS  fcllcv»s: 

INPUTS:  PAYLOAD,  ALTITUDE  OPTION  KEY,  ALTITUDE.  FILM  TYPE  AND 

THICKNESS.  TAPE  TYPE  ANC  LOAC  RATING.  NUMBER  OF  1A?CS,  AND 
CAP  THICKNESS 

OUTPUTS!  NUN-C  IMtNS  IONAL  FILM  AND  TAPt  HEIGHTS.  TAPE  STIFF- 
NESS, AND  FILM  MODULI. 

I******************************************************************* 

ALL  EMFIRICAL  RE L AT  I CNSH I PS  USED  ARE  SUdJECT  TO  REVISION  4HEN 
AODIIIUNAL  TEST  DATA  dECCME  AVAILABLE. 


THIS  SUBROUTINE  ASSUMES.  FCR  LACK  OF  A GCCC  TEMPERATURE  MO OFL , THAT 
THL  BALLOON  MATERIAL  1 _ MPE  RAT  URE  IS  ATMUSPHERIC  TEMPERATURE. 

THIS  ASSUMPTION  MAY  NOT  BE  ADEQUATE. 


HEAL  K.KT.KTAPE, LAMBDA 
REAL  MG.MA.LF.IAM 

INTEGER  CCDEF.CCDET 


ICC  SERIES  - DEFINE  NECESSARY  VALUES 


IF(KEY.EQ.2)  GC  TC  101 

100  H =LUN  ST 

H = F*0  .3040 

CONVERSION  OF  H FROM  FEET  TC  METERS  FCR  SUERT  ATMOS 

CALL  ATMCSIH.TA.PA.RHCA.E. IAM.GXPAN.TOHICH) 

GU  TC  1 C2 

101  P A =CC  NS  T 

CALL  ATMCS2 I FA ,H , TA.RHCA. B.  I AM , GXP AN , T OH  1 G F > 

102  B = d*fc  .243E-02 

CCNVERSI ON  CF  0 FROM  KG/M3  TO  LEM/ FT  3 
T=TA-273 . 1 1> 

CONVERSION  UF  TEMPERATURE  FRCM  DEG(R)  TO  GEG(C) 

LAMBDA  = (P/B  )**( 1 ,/3.  ) 

K = 1 ./((  6.263185)**  (I. /3.  I) 


200  SERIES  - COMPUTE  FILM  PROPERTIES 


IF < CLCEF  .EO  .2  ) GO  TO  202 
201  KFILM  = 5e. Cfc*<  F THICK/12  . ) 

EFILM  = ( 0.02*< T**2 )-l .60*1 T )F41 . 951*1000. 

GC  TL  203 

2C2  mFILM  = 87. CE*<F THICK/12. ) 

LFILM  = 777??  MYLAR  FILM  DATA  IS  NECESSAkY 

203  SIGMA  = mFILM/(K*U*LAM0DA  ) 

CSIGMA  = SIGMA*  < CAP/ FT  Hi CK  ) 

C 

C THIS  PROGRAM  ASSUMES  A POISSON'S  RATIO  UF  D.50 
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.EVEL  21  MAIL  DATE  = 78226  21/29/42. 


fcV  = CfcF  ILM*FTH  ICKA12  . >/ ( B *(  L AMECA**2  ) > 
EC  = EM 
EMC  = EM/2  • 


^ 300  SERIES  - CCMFLTE  TAPE  PROPERTIES 

C 

I F (CODE T.EG.2 > GC  TC  302 

301  * T AP  t = ( E . 1E-CC  MTLW+  C.  CC38*  . 3024*  (.  003*F  THICK  ) 
»TAPE2=*TAPE*.3024*(CAP-ETEICK  ) 

K T APE  = ( TLR*10.  > *(  1 . - ( b . 6 7E  - 03  > * (T-20.  ) I 

G J TO  303 

302  *TAPE  = (2.0E-06 >*TLR*0.0021 
K TAPE  = 26  • * TLH 

C 

C NC  INFORMATION  CN  LC»  TEMPERATURE  EXPECTS  CURRENTLY  AVaILALLE 

C 

303  TSIGMA  = (NT  *WT  APt  >/ ( E *(  L AMBl>A**2  > > 

TS  1GMB=  TSIGMA*VkT  APE  2 / A T A PE 

K.T  = (NT ♦KTAPE  ) / ( B * ( L A M BO  A * * 3 ) ) 

RETURN 

END 


oonOnrtr.noAr'finf'rirn 


-EVEL  21 


A7MJS 


DATE 


78226 


? ! /^9/42 


SUERCLT  INE  ATMCS ( H . T A • P A , R HO A . B • IAM, GXPAN.TOHI GH) 

......  SUBROUTINE  FUR  SOLVING  FOR  THE  VALUES  OF  TEMPtRATURE. 

PRESSURE.  DENSITY,  SPECIFIC  HUUYANCY,  INTEGRATED  AIR 

......  MASS.  ANO  "GAS  EXPANSION"  FOR  ANY  GIVEN  ALTITUDE  BELOW 

......  61000.  METERS  <2001.11  . FEET).  ALL  E0UAT1UNS  HAVt  DL.C.N 

......  DERIVED  ACCORDING  TC  THE  U.S.  STANCARL  ATMOSPHERE.  1962. 

ALL  UNITS  ARE  IN  THE  SI  SYSTEM.  I.E. 

TL MPERATURE < CEG  KELVIN),  FRESS'JRE  (MO)  « CtNS  I T Y<  KG/M  3 ) , 

SPEC  IF  I C BOUYANC YCKG/M3)  INTEGRATED  AIR  MASS(KG/M2>. 

GAS  EXPANS  IUN<C  IMENS IONLESS  ) . LIFTS  WERE  BASED  ON 

GRADE  A HELIUM  AND  FJRF  HYCRC-CEN. 

FIR  SPECIFIC  eOUYANCY,  TC  CONVERT  FROM  Ko/M3  TO  LH/FT3. 

MULTIPLY  tit  0.06  24  3 

H IS  INPUT  IN  METERS 

REAL  MG. MA.LP.  IAM 

T0H1GH  =1 
MC=4  . CC26 
MA-26  .9644 
R=E. 31432EC3 
GY=9 .806b5 
RH0S1 F = 1 .2250 
I GA  5=1 

IF  HYDROGEN  GAS  IS  USEC  INSTEAD  OF  HELIUM, IGAS  MUST  BE  CHANGED  TO 
SOMETHING  CTHER  THAN  1 

IF < H .LT .olOOO  .0  ) GO  TC  10 
WRITE  <6.155)  H 

155  F JRMAT < '0 • . 5X. • ♦***•  ALTITUDE  H UUTSIOE  RANGE  OF  SUBROUTINE  ATMOS. 

1 H =*,E14.7.»  METERS.  *****♦•) 

I OH IGH=- 1 . 

GO  TC  28 

10  IF  <H.LT.52CCC.0)G0  TC  12 
HE=  52  CO  C .C 
T E =2  7 0.65 
LP=-C.0C2C 
P0=O .590 
RH0t)  = 0. 0007594 
GJ  TO  22 

12  1F<F.LT  .47000  0 ) GO  TO  14 
HB  =4  7 C 0 C • C 
T E=  2 7 0.65 
LP=0.0 
Pc3  = l . 109 
RFUH=0  .00  142  75 
GC  TC  22 

14  1F< H.LT. 32000. 0)  GO  TO  16 

HE  = 32 00 0 .0 
Tfl=2?€. 65 
L P=  0 . C 02  8 
P tj  = 8 • 6d  0 
RHOH=  C.  Cl  3225 
GO  ro  22 

16  1F<H.LT.  20000.0  )GC  TC  18 
HB =20000. C 
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TB=216.65 
LP  = 0 .00  10 
Pa=54 . 745 
kt-OB=  C .086035 
GO  TO  22 

18  IF  ( H .L T . 1 1 CCC. 0) GO  TC  20 
He=  1 1 000  .0 
T 3=2 1 6.65 
LP=0.  C 
PE  = 22  6 . 320 
RHOB=  C. 36352 
GO  TO  22 
20  Hb=0.0 

Tb=2fi  6.  15 
LP=-0  .0065 
P 0 = 1 0 1 3 • 2 50 
«HOB=  1. 2250 
22  T A=TB+LF*(H-(-E  1 

I FILP.EG. C> GC  TC  24 

P A=PB/( LP*(  h-HB  J/TB-M . C) **{ G V*MA /L P/R) 

GC  TL  26 

24  PA=PB/E >P< G Y* VA/P/ TB* < H-Ht > > 

26  «HUA=  PA*MA/k/TA *100.0 
I F ( IG AS . NE . 1 ) VG-2.J1 55 

b =RHG  A * ( 1.0-M3/MA) 

1AV=10.0*PA*(1  .03  75  1-  0 .00527*ALUG10< PA  ) ) 
GXPAN=HHGSTP/hHOA 
2b  METURN 
END 
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LEVEL 


? I Z29/42 


21 


ATMOS  2 


DATE  = 78226 


SUURCLTINE  ATMCS2 (PA  ,h  .T A .RHCA  .6  . I A V , GX P AN. T CH l GH ) 

..  SUBROUTINE  FOR  SOLVING  FOB  THE  VALUES  OF  TEMPERATURE. 

..  ALTITUOE.  DENSITY.  SPECIFIC  ECJYANCY.  INTEGRATED  AIK 
..  MASS.  AND  " GA  S EXPANSION"  FOR  ANY  GIVtN  ALT  I TLDE  QELOVk 
..  61000.  METERS  (200131.  FEET).  ALL  EQUAT 1UNS  HAVc  BEEN 
..  DERIVED  ACCORDING  TO  THE  U.S.  STANDARD  ATMOSPHERE,  1962 
.. ALL  UNITS  ARE  IN  THE  SI  SYSTEM.  I .t  . 

..  TEMPERATURE (DEG  KELVIN).  PRESSURE (M3),  CE NS  I T Y ( K G/M 3 ) , 
..  SPECIFIC  UUUYANCYt  KG/M3)  INTEGRATED  AIR  MASS(KG/M2>. 

..  GAS  EXPANSICNCCI MENSICNLESS  ). LIFTS  HERE  3ASED  ON 

..  GRADE  A HELIUM  AND  PUKE  HYDROGEN. 


FCP  SPECIFIC  HCUYANCY. 
MULTIPLY  t>Y  0.C6243 


TO  CONVERT  FRUM  K G/M 3 TU  LB/FT3, 


PA  IS  INPUT  IN  Me 

REAL  MG.MA.LP.  IAM 
T UH IGH= 1 
MC  = 4 .0026 
M A =2  S . 9664 
R-8.31432E03 
GY  =9.80665 
KHOSTP-1 . 22  S C 
I C A S=  1 

IFCPA  .LT .226 .32  ) GC  TC  10 

Pb=1013.25 

FG=0  .0 

T0=2«a. 15 

LP=-0 .0065 

PHCO= 1 .225 

Go  TO  22 

10  IF (PA  .LT  .54 .749  ) GQ  TO  12 
PE  =2  2 6 .32 
HH=1 1 CCC. 

T B = 21 6 .65 
L F =0  • 0 

«HOB=  C. 36392 
GL  TC  22 

12  IF(PA.LT.8.680)  GC  TC  14 
P E=  5 4 .749 
H6=20000  • 

TU  = 2 1 6.  65 
LP=0 .00  1 
KHCB=0. 088035 
GJ  TO  22 

14  1 F ( PA  .LT  . 1 . 1 09  ) GO  TO  lo 

P8  = 8. 6ec 

HE=  32C0C  . 

T 0=228.65 
LP  = 0.  C 02  e 
RF0B=0  .01  3225 
GC  TC  22 

lfc  IF  ( PA  .LT.  C.  590)  GC  TO  18 
P E = 1 .109 
HU=47C0C. 

T P=  2 7 0.65 
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XV  EL 

la 

20 

155 

22 

24 

26 

28 


21  A T MGS2  OAT  E = 78226 

LP=O.C 

RhCfl=0. 0014275 
GO  TO  22 

1F(PA  .LT  .0 .1828  ) GO  TO  20 

pa-o. sy 

Hd  = 52  coc. 

T E =27 0 .65 
LP=— C.0C2 
R PQ B=  0 .G0075S4 
GO  TC  22 
MRITEC6.1S5)  PA 
WR1TEC6. 155 ) PA 

FORMAT! *0*  .5X»* *****  ATMOSPMFRIC  PRESSURE  04  OJTSIDE  RANGE 
10UTINE  ATMOS2.  PA  =*,E14.7.'  MH . **•**•) 

TCP  IGF=-l  . 

GO  TO  28 

1F(LP  .EQ.0.0  ) GO  TO  24 

H=PU*TB/LP* ( (P6/FA)**(LP*R/GY/MA )- 1 .0 ) 

GO  TO  26 

P=ALCC!  FB/PA ) *R*T E/GY/M A*FR 
TA=TB*LF*  1H-H8 ) 

RHOA=PA*MA/R/TA*100.C 
1F(  IGAS.NE.l  ) MG  = 2 .0159 
d=RHUA*! 1 . C-MG/MA) 

I AM= 1 0 .0*PA*(  1 .0 375 1 -0. 00527* ALOG  IOC  PA)  ) 

GXFAN=RHCSTF/RHCA 

RETURN 

ENC 
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FJBR 
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-EVEL  21 


QO  VNCV 


l)A  TE 


78226 


? 1/29/42 


I 


SUBROUTINE  atVKCY  (B  .PA  . I GAS  .TOHI  GH  ) 


SUBROUTINE  TC  SCLVE  FCR  ATMOSPHERIC  PRESSURE  (PA)  GIVEN  T PC 
CORRESPONDING  GAS  BOUYANCY  <B)  AT  ANY  ALT1TUDF  OELCW  61.000 
METERS  (200.131  FEET).  ALL  EQUATIONS  HAVE  BEeN  DERIVED  ACCORDING 
TO  THE  U.S.  STANDARD  ATMOSPHERE.  1962. 

ALL  UNITS  AWE  IN  THE  SI  SYSTEMS 


SPECIFIC  UUUYANCY  (B)  IS  THE  INPUT  IN  KG/  M3 
ATMOSPHERIC  PRESSURE  (PA)  IS  THE  OUTPUT  IN  MB 


CALCULATIONS  ARE  BASED  ON  PURE  HYDROGEN  AND  GRADE  A HELIUM. 
FCK  HYDROGEN.  INPUT  IGAS  AS  0. 

FOR  HELIUM,  INPUT  IGAS  AS  1. 


TCFIGF  WILL  EE  SET  TO  -I  IF  BOUYANCY  IS  OUTSIDE  THE  PARAMETERS 
OF  THE  STANDARD  ATMCSPHEWE. 


REAL  MW 
T CH  IGH=1 
M W = 4 . C026 
B 1 = 1 .0557 
02=0.31 355 

0 3 = 0. C75e58 
B4=0 .01 1396 
35=1 .2308E-3 
tJ6  = fi.  E44  4E-4 
07=2.36356-4 

1 F ( IGAS. EC. 1 ) GC  TC  90 
MW=  2.01 59 

ftl =1.1397 
t)  2=0. 33e53 
83=0 .081899 
04=0.012303 
05  = 1 . 32eeE-3 
Bo  = 7 .0650E-4 
37=2. £51 7E-4 
90  IF(B.LT.B1  )GU  TO  1 
• HITE  (6  .101  ) E 

101  FORMA T(  • O'  • ***********  VALUE  OF  E IS  OUTSIDE  LIMITS  OF  ST4NCAP0  AT 
1MOSPHERE.  B = *.E12.5.'  (KG/M3)  ****♦**¥**•) 

TCHIGH=-1 • 

GO  Tf  40 

1 IF (b.LE .82 )G0  TO  2 
T0  = 2dfc.  1 5 

XLP=  — C • C 065 
PE  = 1 013  .250 
GO  TO  9 

2 1 F ( B .LE  .83 >GU  TO  3 
TB=21 6.65 
XLP=O.C 
PB=226.3 20 

GG  TU  9 

3 IFCO.LE  .84  )GU  TO  4 
TB=21 0.65 


-EVEL  21  3UVNCV  DATE  = 782  2b  21/29/42 

XLP=0 .0010 
P8=S4.745 
GO  TO  9 

4 l FIB. LE.Bb) GC  TC  5 
TB=22€. 65 
XLP=0  .0028 

P8-=  8 . eec 

GO  TO  9 

b IF (B.Lfc .86 )G0  TO  6 
T8=27C. 05 
X LP=  0 .0 
P8  = l . 109 
GJ  TO  9 

o IF<B.Lfc.87)GC  TO  7 


T8=2?C.65 
XLP=- C • C 020 
Pe=0.590 
GO  TO  9 

7 WR  lTfc  <6. 101  1 8 
T OH I Gh  = - 1 . 

GO  TO  40 
9 ITTERsO 
PNE»=F8 
10  PULD=PNEW 

ITTER=ITTER*l 

FPA=POLD-<  2. e7C531*  fi*  TQ/  ( 1 . -Vw/28  .9644)1*11  .♦<<  Pti/ P JL  D ) * * < 29  .271  J* 

ixlp  >-i . n 

GPA  = l.*(84.024l*E*TE*XLP*FE/((l  .-MW/ 28 .9044  ) *PULD**2  1 ) * ( P8/P0L0 > ** 
1(  29.271  34XLP-1.  > 

PNfcW  = FCLD-F  F A/GP A 

IF(A6S<  <FNEW-FCLC)/FNEw>  .LE.  .00001  )GC  TC  20 
l F ( ITTER  .GE  .5C  )GO  TO  2C 
GC  TO  10 

20  OELPA  =PNE  W-PCLO 

WR  IT  E ( O t 1 00  ) OELPA 

100  FCGMATI '0 • . 10X  , 'SCLUT  I CN  CIC  NCT  CONVERGE  AFTER  bO  I T TER A T I O NS . • 

1/oX  . 'DIFFERENCE  BETWEEN  NEW  ANC  OLD  FRESSUWt  CN  LAST  ITTER AT  ION  WA 
2S  ' • 2X  . E 1 0 .3  > 

T CH IGF-  — 1 . 

20  PA=PNEW 
40  RETURN 
END 
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lCVEL  21 


GLNGTH 


DATE  = 76226 


21/29/4? 


SUBROUTINE  GLNG  TH( GP  * SUMO A . 1 NO  1C) 


SUBRCUTJNL  TC  SOLVE  FOR  THE  GCRE  LENGTH  RELATIONS  OF  ZERG-PRE SSU& 
NATURAL -St-APE.  FLAT-TOP  BALLOONS.  TO  OBTAIN  GORE  LENG TH/L  A MBl) A 
(SLMOA)  FRCV  GRCSS/FAVLCAO  (GP)  IMT1AL1/C  INOIC  TO  1.  TO  OBTAIN 
THE  REVERSE  OF  THE  ABOVE.  1NT1AL1ZE  INOIC  TO  ZERO. 


THE  FOLLOWING  IS  A CURVE-FIT  CF  FIGURE  5 IN  NC AP-TN/ 1 A-99 • SECT  10 

1F(  IKC1C.EQ.0  ) GO  TO  60 
IF(GP  .GE  . 1.  )GC  TC  10 
AO  WR IT  E (o  » 1 00  1 GP 

100  F CHMAT( • 0* •• **********  VALUE  OF  GROSS/ PAYLOAD  OUTSIDE  RANGE  OF  SUB 
1R  OUT  I NE  GLNGTH.  GROSS/PA YLOAD(GP)  a • .E12  .5  .* ***********  ) 

1 NC I C *-l 
GO  TO  99 

10  1F(GP .GE .1 .1458 )GO  TO  11 
A>1 .99442 

B a C • 2 74  64  9 
GO  TO  50 

11  I F (GP .GE. 1 • 3205 )GC  TC  12 
A ai .95267 

B = 0 .260545 
GL  TC  50 

12  IF( CP  .GE .1 .5306  >gO  TO  13 
A a 1 .969503 
6=0.266635 

GO  TO  50 

13  1 F (uR.GE. 1 . 7630 >GC  TO  14 
As 1.96507 

6=0 .291661 
GO  TO  5 C 

14  IF(GP  .GE .2 .0697  )GU  TO  15 
As  1 .97653 

a=c.2975ei 
GC  TC  50 

15  1F(GP.GE.2.4566)GC  TC  16 
Aa  i .971  75 

B =0 .3022  39 
GC  TO  5 C 

16  1 F ( GP  .GE .2.9C35  ) GO  TO  17 
A= 1 .96326 

6=0.307  C2  3 
GO  TU  50 

1 7 I F (GP .GE.3.4J80 )GC  TC  16 
A = 1 .95482 
8=0 .3 1 1076 
GC  TC  50 

16  I F ( GP  .GE .4.0776  )GO  TO  19 
A= 1 .94521 
B=C. 315067 
GO  TO  50 

19  I F (GP .GE.4.8360 ) GC  TC  >0 
A = l .93701 

8=0.318071 
GC  TO  50 

20  1 F ( GP  .GE  .£  . 7363  )GO  TU  21 


I 
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A=1 .92864 
8=0.320819 
GU  TO  50 

21  I F ( CP  .GE.6.70981GO  TO  22 
A = 1 .92034 

8=0.323266 
GO  TO  50 

22  IF<GP.GE.8.0155)GC  TC  23 
A=l. 91325 

0=0 .325220 
GC  TO  5C 

23  1F( GP  .GT .10  . >GU  TO  40 
A = 1 .9  066  7 
0=0.326675 

50  SLMOA=A*GP**U 
GC  TO  99 

60  IF  ( SLMOA  «G£  • 1. 99442)  GO  TO  70 
91  Wfc ITfc <6 .200 > SLMOA 

200  FORMAT! • C* .* **********  VALUE  CF  GORE  LE NOT F/L AM80 A OUTSICE  *» AN  Gt  O 
IF  SUBROUTINE  GLNGTH.  SLMOA  = ••£  12.  E»* **********'  ) 

INLilC=-l 
GO  TO  99 

70  1F1SLMOA  .GE  .2 .07044 )GO  TO  71 
A = 1 .99442 

8 = 0.2  74  649 
GO  TU  95 

71  IF(SLMDA.GE.2.I5453  ) GC  TC  72 
A =1.9 92 67 

B=0 .280545 
GC  TC  95 

72  IF C SLMOA .GE  .2 .24767  )GO  TO  73 
A = 1 .989503 

8=0.266635 
GU  TO  95 

73  I F (SLMOA. GE. 2. 35038  ) GC  TC  74 
A=  1 .96507 

8=0.291881 
GO  TO  95 

74  IF ( SLMOA .GE  .2  .46373  )GU  TO  75 
A = l .97853 

8=0.297561 
GC  TC  95 

75  I F ( SLMOA. Gc  • 2 .58781  >GC  TC  76 
A= 1 .9  71 75 

8=0.302239 
GC  TC  95 

76  1FTSLMOA  ,GE  .2.723391GO  TC  77 
A = 1 .96328 

8=0.307023 
GO  TO  95 

77  I F < SL  MO A .GE. 2. 87037  ) GC  TC  78 
A = 1.95482 

B=0 .311076 
GC  TO  95 

78  1 F ( SL  MOA .GE  .3 .02694  )GO  TO  79 
A = 1 .94521 

8=0.315067 
Go  TU  95 
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79  IFtSLMOA.GE 
A = 1 .93701 
8=0.318071 
GC  TC  95 

80  IFCSLFDA.GE 
As  1 .92864 

6 = 0.  320819 
GO  TO  95 

81  IF(SLMOA.GE 
A =1  .92034 
B=0.3232ee 
GC  TC  95 

82  IF(SLNOA.GE 
A= 1 .9 1325 
8=0.325220 
GO  TO  95 

83  l F (SL  MOA  »GT 
A = 1 .906  6 7 
8=0. 326675 

95  GP=(SLMDA/A 

99  RETURN 
END 


3.19819) GO  TO  80 

3.37780  )GC  TC  81 

3 .56703 IGO  TO  82 

3.76484 )GC  TC  83 

4.0  )GO  TO  91 
* * < 1 ./E> 


APPENDIX  B 


Balloon  Design  Program  - Input 

The  input  data  for  the  program  consists  of  three  cards  per  balloon  to 
be  designed.  In  order  to  stop  the  program  without  generating  an  error  state- 
ment a similar  set  of  three  cards  with  no  data  is  required. 


Card  1 . 

Format 

(F10.0,  110, 

fio.o,  no,  fio.4,  no,  fio.o,  no) 

Col umns 

Variable 

Description 

1-10 

P 

Design  payload  in  pounds 

20 

KEY 

Altitude  option:  1 -Altitude  in  feet 

2-Altitude  in  millibars 

21-30 

CONST 

Design  altitude  in  feet  or  millibars 

40 

CODEF 

Film  type:  1 -Polyethelene;  2-Mylar 

41-50 

FTHICK 

Film  thickness  in  inches 

60 

CODET 

Load  tape  type:  1 -Polyester;  2-Kevlar 

61-70 

TLR 

Tape  load  rating  in  pounds 

71-80 

NT 

Number  of  load  tapes 

Card  2. 

Format 

(F10.2,  F10. 

2,  F10.2,  F10.2,  110,  F10.2,  F10.2,  F10.4) 

Columns 

Variable 

Description 

1-10 

TL 

Top  load  in  pounds,  (+)  up,  (-)  down 

11-20 

TAUO 

Stress  constant,  0.  for  natural  shape 

21-30 

TAUI 

Stress  constant,  0.  for  natural  shape 

31-40 

ALPHA 

Superpressure  in  feet 

41-50 

N 

Print  increment,  0 for  all  points 

51-60 

DSO 

Nondimensional  gore  increment  (ds/x) 

61-70 

CSTART 

Nondimensional  gore  position  of  edge  of 
oap  (Scap/x) 

71-80 

CAP 

Film  thickness  including  cap  in  inches 
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Card  3.  Format  (3110,  2F10.0) 


Columns 

Variable 

Description 

10 

KEY2 

Output  control  (0-print  only,  1 -punch  deck 

& print,  2-disk  file  generated) 

11-20 

MPT 

Number  of  points  desired  in  load-altitude 

curve 

30 

IGAS 

Identify  lifting  gas,  1 -Helium 

31-40 

PMIN 

Minimum  Recommended  payload  in  pounds 

41-50 

PMAX 

Maximum  Recommended  payload  in  pounds 

The  output  of  this  program  consists  of  a),  the  nondimensional  shape  para- 
meters normalized  with  respect  to  A;  b).  the  dimensional  manufacturers  table 
coordinates;  c)  100  equally  spaced  points  needed  for  input  into  the  analysis  program; 
and  d).  the  load-altitude  curve.  A sample  of  the  printed  output  is  as  follows: 


OOOOOO  300000 OOOOOOOOOOOOO OOOOOOOOOOO 0000  00  090000 0000  90000 
oooooaooooooooooooooooooooooooooooooooooooooooooooooooeoo 
OOOC  OOOOOO OO OOOOOO OOOOOOOOOOO OOOOOO OOOOOOOOOOOOOOOOOOODOO 
ooo  oo  ooooe  oo  oooooooooooooo ooooooooooo oooooooooooooooooooo 


53553«55i;5553535555555553553555555313555555555555S555555 


0*00>0000^--^CVNN.>JI») J>4tA«Ain®  $ O <> « O <0  O <0  K KK 

ooo ooo oooooo oooooooooooooo©©© ©oooooo OOOOOO OOOOOO OOO OOOOOO 
• • • • • • 
oooooooooooooo  ooooooooooo ooooooooooo ooo ooo oooo ooooooooooo 


iiiiiiiiilisiiiiiliiiiiiiiilsisiiiiliiliiissiiiisisiiilii 


2SiSS;?S*3l335?;5=?SS3SS5SS|S|aS2|||5S8S«SS«S|?i«|«|82jg2 

OMfOCON*  ASON«OSOMrtn<)AO**Nn  l>8  8NB8®*i8N'0tflrt»»^N8OdMKM ®OOi«<©K®^© 
• ♦••••••••••••••••••••••••••••••••••••^  • • ••••••••  •••••••• 


- M 10  ♦ rf)  0 N * 9 O m»  2 ♦ <2  o N • g O £ Jj  JJ  <*  jg  ^ {J  • !5SS9!3S<1^3SS?;33:5235?S3aSSSSK 

J*3j*5o5ooOOO  J JJJ  JvJJOOJ  JJO 


<0^«onNMB9>nN**B»N®9tBN8O«  O-rflOlOK  — lOO'OOO^OMOO^ONAO’OK-AOSfcr^S 

• • • • •••••••  • • • 


Mi  -«“»*-»•— ••2=2!2S22t:22Ss:a3SSSS;5RSS»SiJ5SiS?R*>t?5058S8^5S3{S3l3*S*5J 


000000000009900003000d0000000900090000e0000 
99OO9OOOOO9O9O9O93O99OO9O9O999O999J#O099O9e 
390000000 >090000 0300 00000009000000 OOO 9 99000 
• ••••••••••••••••••••••••••••••••••*•*••••• 

90000000099900000 30  000 0000000900  00000000090 


SSW33SSJ$M8KB8«ttK3.M2:2*S*S?S38m8;IM 

335j3s»3 


6 000  00  0000  00  0000000000  90000 


0009000000000000 

0000000000000099909000900000090009090000009 


V)<0O999OO9  0 •••••• 

liiliiiSllIsillllillililiilislillslIlillli 


2!*5?!cl51:55;:sS?!s555!53i*?53e,5:5s?25s??I 

9 >«  SKf*9  Oi>*l9N«*0»OP>iA#  9N09N9AflN09^«9<9**0«^i99  H — O 


SSS32?2525SS8:22Sfi8C22§;222S85SS?;855S#S8*o 


535 5533533 5 335 35 55 5555 555555o5o 5oo 
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LUAD  - AlIIIUDL  DATA 
<iKU5S  AlRbOWN 

WL  l OH  1 - t Kt> ) AbTlTUDb  - (AM) 


0 

.2386b 13E 

04 

0 

3045546E 

02 

0 

• 1 Cl  L 

04 

0 

304 2964  l 

02 

0 

•2409509b 

04 

0 

304  04  0 2L 

02 

0 

• 2420b69b 

0« 

0 

303 7b 1 9b 

02 

o 

• 24  322076 

04 

0 

3036216b 

02 

0 

• 24  4.ai»50E 

0* 

0 

3032632b 

02 

0 

.24t>4404L 

04 

0 

3030070b 

02 

0 

• A4bb<;bi;l: 

04 

0 

30776096 

02 

0 

.?47  7b0?b 

04 

0 

5024  9t>biL 

02 

0 

• i 4 bb9  50E 

04 

0 

30224c6t 

Oc 

0 

•250079bb 

04 

0 

301 9906L 

02 

0 

.25  1 1 647b 

04 

0 

30 1 73b 4 L 

02 

0 

,ia7i:v95L 

04 

0 

301 4bb6t 

02 

0 

34  5E 

04 

0 

J012406L 

02 

0 

•2545693b 

04 

0 

300  9926b 

02 

0 

.2  55  7 04  IE. 

04 

0 

300  74  4 6b 

02 

0 

04 

0 

joooooat 

02 

0 

• 26  79738b 

04 

0 

3002560b 

02 

0 

.2691 086C 

04 

0 

30001 33b 

02 

0 

• 2bO«^4  36b 

04 

0 

299 7694 L 

02 

0 

• ?t>  1 3/b  4b 

04 

0 

29952976 

02 

0 

•2b261  33t 

04 

0 

2 99..tltt0b 

02 

0 

• 2t»364t-lE 

04 

0 

2990 60 4 t 

02 

0 

.764Tb  31P. 

04 

0 

2988129b 

02 

0 

.26591 79b 

04 

0 

2985763b 

02 

0 

• 26  7 0 52  7b 

04 

0 

2983397b 

02 

0 

•26blb75b 

04 

0 

298104  lb 

02 

0 

•2o93?24fc 

04 

0 

t97H706b 

02 

0 

• 2 704674b 

04 

o 

2976372b 

02 

0 

.271 5922t 

04 

0 

2974069b 

02 

l) 

• 2 72  72  70L 

04 

0 

2971 74 4 t 

02 

0 

.273801 6b 

04 

0 

2969460b 

02 

0 

•2749967b 

04 

0 

2967167b 

02 

0 

.27612  lot 

04 

0 

29o48b6b 

02 

o 

•2772665L 

04 

0 

296201  lb 

02 

0 

• 27640  1 ol. 

04 

0 

2960369b 

02 

0 

. 2796  261 L 

04 

0 

2958107b 

02 

0 

•280b709b 

04 

0 

2966836b 

02 

0 

• 2ol  6059b 

04 

0 

2963561b 

02 

0 

• 2829*.  Q8L 

04 

0 

2961309b 

02 

o 

.2fc  4 0 756b 

04 

0 

29490596 

02 

APPENDIX  C 


Balloon  Analysis  Program 
BALAN 


//a  IONS 

C**  * **  ****  »***  ***********  ******************************  ************* 

C * b A l A N — 

C tlAcLUON  ANALYSIS  PHUOWArt 

C* ***************  *************************************************** 


***< 


***« 


C This  PKOCiRAM  IS  i NT CNDLO  TO  CuMPUlt  1HL  SHAPE  AND  SIMESS 

V_  DlSTKluUT  IONS  OP  A hALLilUN  AT  ANY  AcTITUDL  •HR kb  Trtb  R ILM 

C MAY  bt  CONSIuLRcO  TO  LOOP  rttTWtfc'N  T APh S « THc  MANUFACTURED 

C kAUIuS  AND  Wl  lGHT  DISTRIBUTIONS  ARP  iTPUUlHtD  INPUTS.  THt 

C Dl  SION  SHAPt  PAHAMhTLRS  ARE  NE  E D t-  D 10  SlART  THC  ITERATION 

C PRwCPES.  THE  OUT  HOT  CONSISTS  OP  THE  OhFOKMlD  SHAPE  PARA- 

C MLlO'ia*  MbklDlUNAL  SlktSS.  C I PCUHI'L  FlfcM  I AL  STRESS.  ShC  AW 

C STEirSS.  PRINCIPAL  STRESSES*  AND  1INS1UN  HtLU  PATTERNS 

v_  AHtnc.  APPL1CABLL. 

C ACL  QUANTITIES  ARb  NUND  l Mt.NGl  UN  AL  1 ZFD  «11H  R^SPt-CI  TU 

C SMAA,  Trtc  JNuLFuRMLD  GONE  LENGTH,  AND  PD,  THl  Dt  ISON 

C PAYLOAD. 

C*********  *********  ************************************************* 

IMPLlcll  Rt  AL  *1.  ( A-H  »U—  /.  I 

DIME  NS  1 ON  A(  oOO  ) ,C<  JbOl  >,b  ( lb01).X(  1 r,01  ) ,k 0(300  ) . WT 
1 ( J00 ) .GP( 30  0 ) .WIAPl  (oOO  >.SmP(  300  ) .CLP(JOO) 

wp.al*e  nu.koar .lamhDa 

AT  AN  I / J -DAT  AN ( t ) 

ARSlNlZ  >=DARS1N(2  > 

AWCOSI l ) —DA KCOt ( t ) 

Si  N C i ) =t)S  l N ( 4 > 

LJSC i ) =DCUS(Z  ) 

1 AN(  Z )-D7AN( Z ) 

kE  AD(S.nOO)  Pu.PIUP.WliaL.SMAX  .LAMi.o  A 
kL  Aj  l j.tOl  ) .lUtXHAH  ,KiiAr 
He  AO  ( s»  * b 00  I < MmX  • Au  AW  ,ubAk 

XLAii(j|K02  ) < GPlK  > • PO  (K  ).WT(K  ),«TAPL(K)  .K-l.XMAX) 

b*AX -GPlKMAX ) 

Si,C=  SMAX/LAMntlA 
IMAX=t>*KMAX*l 
1>U  10  K-|  .KMAX 

L-0*(A-1 1*1 

kPADlS.nOJ)  I X I L* 1 ) . 1-2. o> 

A(  c* 1 ) -0 .0 
eON I INUL 

Hi  -4.*A1A  4(  1.0G0) 

PuN— PI /N 
Al 1 >=1  ,-PTUP 
ttoAH  =t*t»AR*SUL  **  3 
Kl»Ak=kHAR/L  ./P  I 
c0AK=rbAx  *SUL 
VTOP-Xl  1 MAX  ) 

DO  JO  K — 1 . X MA  A 
C=6*|K-| )*1 
X(  L*2)-X  (C*«.  >/SUC 
A(L**>)  — X(L*31^S  ol 

XI  L*o)  - I V T JP-  X I L*t>  > )/SUL**j*PI 
GP(K  )=oPIO/SMAX 
*0|lO-KOCA.)*N/l<r./la  I / SM  AX 
• 1 IK  )=WT(KI*SUL 
w T APL  I K ) =■  •!  AP  (_  ( K » *SDL 
UIiyT  iRUt 

Mkl I t<  O.W03  X K »GP|K  > . kGIK) .NT  IK  1 , • 1 APE C K ) *K— 1 » KMAX ) 

SPUN— S 1 N ( PUN ) 

CPUN  —C.OS  I PuN  ) 


**** 


1 0 


I 

•! 


71 


* "V*  * 


I 


: 


4 2 
4 3 
44 
4 b 
4 fa 
47 
4tt 
4'* 
bO 
fa  l 
b2 
fa  3 
b4 
:>b 
b6 
fa  7 
bfa 
b9 

00 

01 

02 
fa  3 
t>4 
fab 
06 
07 
fad 
ovl 
70 
/I 

72 

73 

74 


7b 
7 fa 

n 

m 

7-4 

oU 

ol 

fa  3 
34 
fab 
fa  fa 

K7 

fab 

fay 

vO 

PI 

'/3 

•A 


P4 
ifa 
•lb 
»/ 
fa, 3 


fa2PON=  oPON*  hPON 
35= 1 ./ C KMAX-1  1 
2PA=AbAR 

00  3 MNN=l.b 
100  .IU  200  K=1.KMAX 

1MK  ,tU.bO)Li>faP=4.6»EUAP 

tX=ubAk/( 1 .-NU**2 ) 

uC-KM 

hMC=EC*NU 

U = LM*l  U-tMC*42 

L>M=hCVO 

DC=LM/U 

DMC  = - L iAC/U 

L^o*  < K- 1 ) + l 

bTH=blNCX<L+fa) ) 

C I H-CUS  1 X(  L + t-  ) ) 

^fa=SJN(X<L+l) ) 

Ui=CUb(X(L+l)  ) 

Al_FA=AHb  1 N(  bPUN*C1H  ) 

U-17PA.LT.0.)  AlFA— — AL  fa  A 
1FCCTH.L1.0.  > AlFA=-ALFA 
AuL  I =AuF  AIX  (L  ♦ 1 ) 

PH  1 = 2 • * Afat£T 
bbPiJN=falN(  AHL  T 1 

FlM=kU(K  )*UC/DMt+lt)C*UM/l)MC-0«C)*KfcAk 
IF(K.hU.  1 ) !HP  = (Xt  l2>-XCO>l/OS 
Ifa(K.cU.KMAX) ThP=<  X CL  +b )— X CL— 11  I/O b 
IF (K.LU. 1 ) OU  TO  102 
IF (K.cU.KMAXl  OO  TO  102 
THP= 1 X ( L ♦ 11 1-X1L-1 1 1/.../D5 
102  CONTlNOc 

tfaA  = XU  + J)  t-AOAk 

5H=  1 X(L  + <t  l-KHAR*DrtC*AHt  I /PON*  (HOC  K 1 *ObAR * 2PA+ STH*  < *T  C K ) 
1 -*T  Apr  ( M/POU/2  . 1 ) ) /( ROCK  1 +KB AH* DM  ♦ KBAR  *DKC*AB£ T /PON* 

2 HU (K  >* I HP  1 
t.F=X  CL34  ) -HOC  K 1 ♦Si4 
IPlcF.Lf.O. )bM=X|L+4)/kUCK) 
lFlrF.LT.O.  ) t F = 0 * 

FTb=bPON*CU*C lM+bo* CCPUN*C TH4CTH-1 • +bTH*bTHl 
bC  = A Hr.  I /I'JN*  C HO  CK  ) * faM  * t HP  ♦ROCK  )*3BAk 

l*ZPA+blH*(*TCKl-WTAPCCKl/PDN/2.  1 1 

bMPCKl=SH 

A ( K ) =( X l L + 4 1-kUlK 1*SM !*.;•*  PON 

lOb  CONT1NUL 

IF  CbC.L  T. 0.  lfaC=0. 
faLPC K ) = bC 

1 XP1 =1  . +DM*SM ♦OMC  *bC 
L CP 1 =1  . ♦DMC*bM*OC*faC 

F 1d  = SPGN*Co*SH*C1H*  CCPUM-i ■ 1 
F2b=F TO/PON 

IFCLF.t-O.O.  lot  l=AHblNC  — POU*  ST  i-i*K  T APhCK)/bC) 

IFCtl  . tu.O.  1 fa  T U F F = f AN  1 fa  t 1 ) 

IF  (Lf  .1  0.0.  loU  TO  104 

HAT  1 tJ=HU  ( K ) *;>M/  c X Cl  *4  l~WOC  K ) * bM  1 

ToP=bUAR*ZPA*X(  L+21 *hMPl*bPUN4PUN*  ( *TCK  )-WTAPfc  CM 
1 /PuN/2.  >*hlH-faU*SPUN*ClH*Cb-HATIU*WTAPt-  C K»*bTH/2. 
LiUT=CH*(WAf  K.*hC  + U.*C  CPUN*CTH*CTH*bTH*STMl  ) 
if  ( bC.Lo.o*  lb  i uFr=/os»ppp9P<» 

IF  ( bC.  cu.  0.  loo  TO  104 
b roFF= lOP/u Ol 

104  CUNTINUF 


--r  - 


72 


120 

i 2 1 

122 

1.-5 

124 

125 

126 
1 2 7 
1 2tt 
1^S» 
1 30 
151 
1 52 


155 

1 j*i 

1 36 

1 JO 

1 37 

1 _>6 
ijy 
1 4 0 
N1 

1 *4  2 
1 43 
144 
1 45 
1 **b 
14  7 

1 43 
1 4"i 
1 j0 

151 

152 
1 63 
1 j4 
1 65 
1 j6 

157 


ASt(i-<  1-SI  b)(  PHI  > ) *1  SPI.N/UHPI.N  >**«  KlH 
IF  l 2 PA.  t.  0.0. ) SIGN- I . 

IF  l ZPA  .tU.O.  ) GO  TO  101 
5 1 GN  —Z  PA/  jAbb  ( i.  FA  ) 

1 01  AbUfc»N=N*( SI N| 2*PUN)  *5 lGN*ASkG  » 

NN  = 1 

JJ  = 0 

IFIK.Cl.  .1  ) GU  lU  170 
IF ( K .cu.KMAX )gk  lo  loO 
lFtK.FU.2)  GO  TO  loO 
15V  IF  ( K • GC  .2  > J J - 1> 

IF  ( 14  .Gt  « 2 )NN  = 2 
loO  CoNTINUL 

C(L  + 1 >=ATA:M<SlUFF  > 

IF  ( K .Cu.KMaX-  1 >GO  lil  lol 

Cl  L + 2k  J J ) =X  ( L -4  )tNN*Dii*l.TH*LMFl 

tlL  + jkJJ  >=X  (L-3)+Nf«*DS*cTH*kMFI 

C(L+-tkJJ)  —X  <0  — 2 )+Nt4*llS*  l*)100*ClMk:,C*6lH/PlJN*Flb> 

IF  ( J J.  t U.  olGO  TCI  1 6 1 

I F ( K.  »L  o«2  1 GO  1 LI  l'jli 

lol  CtL+5-JJ)=X(L+ll)  kNM*U^/XIL  + 4 )♦  ( W I ( K ) *S  T H ♦ i-jtfAU  * 2 P A*X  ( !_♦  2 )*LMP1 

i*spuk/ijun-sc/2UMH  ») 

C(L  + t-JJ ) - X ( L + 1 2) +NN*L>5  *EMP 1 * X( C + 2 1 *X  CL  + 2 ) * ASUHN* CTH/2 . 
IF(JJ.cG.O)  Gil  TU  200 
IF  < 14  .CU.KMAX-  1 >NN  = 1 
IH4.tU.XMAX-  1 ) JJ=0 
IF  IK  .CC .KMAX- 1 > GO  TO  1 C 1 
GO  TU  cUO 
170  CONTImUC 

Cl  1 >=nnAk*X (?  )-PlUF-WhAL 
Ct  2 ) -A  f AN l 5 TUFF  ) 

C(j)-rVU(  1 ) 

C ( 4 ) -0 . 

C(5)-X( l I/CTM/2./P1 

CC  o)  =X<  12  )4D5/X  f 5)*  <STt1*NI  ( K 1 Ffcitj  Ak*ZPA*X  ( 3 l*EMPl  * iPON/FUN 

1 -SC*F2u) 

Cl  7)=X  l 1 2 ) + L>S*t  MP1*X(  j ) *X(  3 )/7,*Cll1»ASUbN 
GO  IU  ,00 
ISO  CONI 1NOL 

C<  L+ 1 ) -AT  AN l STUFF  ) 

C(L+2)-kU(KMAXl 

ClL  + 3)=X(i_-31-»U6*CTH*tMPl 

c(L  + 4)-ML-i:I  +OS*l  W 1 ( K ) *CTH+SC*6  » h/Vl)N*F  11*  > 

C( L + 5 ) — — P 1/2. 

Pl<=ribAK*ZPA*KUl  I4MAX  >**2  *P  1 
PSUP=PTUP-PK 

IF  (PSOP.GT  . 0.  ICIL+j >=-APCOSl  C PP-P  TlIP  1/2  . /P  l/K  ( L + 4 1 ) 

IF  (PSuP.lT.  0.  )C  < L+6  ) = Ah  COS  ( ( -PP + P 1 UP  1X2 . /P  1/XlL  + A ) » 

C<L+6>=0. 

200  Continue 

tOAb-huAk /4  < o 

» 

51  6MA=  * 50 
DO  210  1=1.1MAX 

F ( 1 ) =C { 1 ) —X ( 1 ) 

X(  I 1=X<  I ) +5  I GMA*F l 1 ) 

PMS=PM5+F ( 1 ) *F( 1 ) 

.10  CONT1MUL 
2 CONTINUE 

OO  1 XL  = 1 »I4  MA  X 

IF  (KL.Cu.  1 > 4NlTc.lO.13)  I4L.XIKI_1.FIKL) 


73 


1 53 


ITE(b.7)  kl  * ( X(  1-m>*Kl-5)  . 1 = 1 .5  ) • (M  I-Fb*KL-5>  • 1 = 1 .5)  • 

l5MP(Kb).5CPlKL> 


159 

1 3 

fUhlrtAT  I15.2FIC.t>) 

lt»0 

7 

r OK  MA  f (IS.UMO.bl 

lol 

1 

CUNT 1NU^ 

1 o2 

CHt  CK=uAb5( F I 1 )/X I 1 >) 

lt>3 

WM3-KMS**  *3 

lb4 

AH  1 T E ( b »*< 00  ) k MS 

105 

MT=OBAR/St)L**3 

1 ob 

VUL=X<  7 ) * SmA X **  3 

lb/ 

PAY- X ( 1 )*PD 

1 Oti 

5 

tON I 1NUL 

1 07 

MK  1 THI  0*703  > PA  Y <LAMbl:A  .BOAR , VOL  • NU  • Act  A R .SMAX. 

1K.UAR  tH  T « N • kMm  A 

170 

V* 

C 

CALL  PK1NT  |KMAX|A»  bMP  * 5CP  .X.WTAPE  . L>S  ) 
r UK  M A T 5 T A TfcMl.  NTS  F UR  INPUT 

. 71 

b 00 

FORM  A 1 I 5l.  1 b.  / ) 

1 73 

•50  1 

FORMAT I 3£ 15*7.115) 

1 73 

302 

FORMAT! 4L lb *7 ) 

17<* 

r.  0 J 

FORMAT  | 5t  15*  / > 

i rt> 

c 

bOb 

FUKMA1  < 1 15.  2t_  15.7) 

V* 

C 

FORMAT  5TArtMLNTS  F UH  OUTPUT 

1 7f» 

L 

-.00 

FORMAT (IX**  WMS  VALUc  = •*fci0.13) 

1 /7 

703 

FURMAT  < 14.4clb.71 

1 75 

->0t 

FUHMAT  < 1H1,  TOX.'ULSIGN  VALUFS"./.  1 OX  . »P  A YL  L’AD=  • .&X.F7.2, 

i • LHS'.IOX,  ‘LAMBDA- • .F4.0.6X.  «dB  AK  - • ,F  7 . 2 . ✓ , 1 OX  • • VULUML=  • . 

• H(ol  ••bX,»NU-,.F3.2,5X.  • ArfAH-*  .F  t *2* /*  1 0 X * • GORE  LtNUlt-'  . 
3F9.2.*  HSIIX.  *KHAK=  * . Fb.  2 *5X*  •hBAH=*  i^btLt/il  OX  • 

4 "NUMBER  Ur  IAPI_S=»  * I5./.10X.  'NOMhEk  l*'  POINTS— ".lA.//) 


1 / 9 

797 

STOP 

1 HO 

fcND 

1 hi 

SUBROUTINE  PkInT  IKMAX.F.-jMP.SLP.X  . WTAPL  .05) 

lb  2 

IMPLICIT  Kt.ALAb  (A-rl.O-2  ) 

1.5 

OIMFNSIUN  F ( 3 00  ) • ScPC 300) .5MPI300 ) .XI  IbOI  ) .wTAPt(300> 

104 

lMAX=o*KMAX+ 1 

1 H5 

47 1Tb I fc .hOl ) 

1 ho 

*k 1 TCIO.hOO  ) ( X ( 1 ) , 1-2. 1MAX ) 

it  7 

AK  I 1 L<  0.900  ) 

lbb 

PI  -4  .*OAT  A 4 I 1 .L»0) 

lhU 

105  K-l.KMAX 

C 

CALCULATE  SHEAR  5THL5S 

1 70 

L 

v_Tli— OCUS  ( X { <>*K  ) ) 

1 71 

IF(K.CU.l)  GO  TO  101 

ll>2 

IF  (K  .F-U.KMAX  ) GU  111  1 02 

1 >3 

TAU-IFl K+l ) — F (A-l )) 74 ./uS-Wl APt <K )*CTH/2. 

174 

oO  1U  103 

1 175’ 

• i 

•1  AO-  (-3.  *-(AI  ♦4.AFI  A*  1 )-F|K*2  ) ) / 4 * ✓ 05— W T APF  IK  ) *C  TH/2  • 

1 7o 

GO  TO  103 

177 

: oi 

IAU-(FIX-2)-4**FIK-i)+3.AF IK) )/4 ■ /OL-WT APL IK ) ACTM7V . 

1 H< 

1 Ob 

3P=5MP(K  )F5CP  Ik  ) 

t. 

CHECK  to  Set  IF  A TLNSIJN  r 1 fc LO  tXI5Tb 

1 77 

I F I j LPI * ) .Lfc • 0.  ) Gl  TU  104 

2 00 

1 AUM AX-uSORl  1 I I 5MPI K )-sCPIK  ) )/2 .00  ) ♦♦<'♦  T AU**2  ) 

201 

•>1  -SP/2  . ♦ 1 AUMAX 

*.02 

52  — 5P/2 .—  T A UMAX 

74 


ANG=L>ATAn(l.GO*TAU/  ( > ) ) /L  . 

lF(bLP(K).GI.  i>MP(K)  ) ANG=ANG*P1  /<_’•  LO 
GU  TO  105 

ANG—  —l)  A I A >K  bMM(K)/l  AU1-IM/2.UO 
bl  -bMR(  K )/JLUb(  ANG)/L)CUb<  ANG) 
b<:  = 0 . 

1 AUM  AX  =artP (x)/2,D0 

m-dltlft.VOl  ) K.F(KI  «SMH(K),  i»CP(K  ) *1  AO.G1  • S/  . T AOMA  X . A N(j 

format  s f ai  t rttNi b fur  outpui 

FORMA  T t cF 10  .4  ) 

F0KMAI(//.4X.»tJLTA,.7X.»H».'»X.*2,.9X.*T».FlX.»THFlA» 

1 »DA, 'VOL • .//) 

f-  lJ  900  FORMATE  1 HI .43X . *MCR .* .OX.  *CIW .* «GX . *MLR .• »4X. «PKl NCI  PAL  bTRfcbSCb* 
1.4X.*MAX*,/.2  9X.*K».bX.*F».7X.*bTKi.SS*.4X.»brHtbS».4X.  •SHtAR  • 

2 » 7 X » • S 1 • *HX»*  bi_  * • 6>X  * • bHkAR  • .bX  . •ANGLE:*.//) 

214  901  FORMAT  ( Jr OX  . 1 A .WFlO.  A ) 

£15  RF  TURN 

21o  lND 


20  J 
2 04 
£05 

c 06  1UA 

2 07 
LOG 

209 

210  l Os 
C 

c 

c 

211  WOO 

212  oO 1 
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APPENDIX  D 

Balloon  Analysis  Program  - Input 

The  input  for  the  analysis  program  consists  of  three  identifying  data 
cards;  a deck  of  manufactured  radius,  film  weight  and  tape  weight;  and  a deck 
containing  the  design  shape  in  nondimensional  units. 


Card  1 . 

Format  (5E15.7) 

Col umns 

Variable 

Description 

1-15 

PD 

Design  Payload  in  pounds 

16-30 

PT0P 

Nondimensional  top  load  (PTOp/PD) 

31-45 

WBAL 

Nondimensional  balloon  weight  (W/PD) 

46-60 

SMAX 

Underformed  gore  length  in  feet 

61-75 

LAMBDA 

Design  A=(PD/bd)**l ./3. 

Card  2. 

Format  (3E15.7,  115) 

1-15 

NU 

Poissons  ratio 

16-30 

KBAR 

Total  tape  stiffness,  NKy/PD 

31-45 

EBAR 

Film  stiffness,  EtX/PD 

46-60 

N 

Number  of  gores 

Card  3. 

Format  (115,  2E15.7) 

1-15 

KMAX 

Number  of  points  to  be  considered 

16-30 

ABAR 

Nondimensional  superpressure,  a/S„„„ 

max 

31-45 

BBAR 

Ratio  of  specific  lift  at  altitude  to 

design  altitude 
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Card  4.  Format  (4E15.7)  KMAX  Required 


Columns 

Variable 

Description 

1-15 

GP(K) 

Gore  position  of  point  K in  feet 

16-30 

R0(K) 

Half  gore  width  at  point  K in  inches 

31-45 

WT(K) 

Weight  increment  between  K and  K+l 
including  tape 

46-60 

WTAPE(K) 

Weight  increment  of  a single  tape 
between  points  K and  K+l 

Card  5. 

Format  (5E15.7) 

KMAX  Required 

1-15 

X(L+2) 

Normalized  design  radius  at  point  K,  r/x 

16-30 

X(L+3) 
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Normalized  design  volume  at  point  K,  V/irX 

Output  consists  of  computed  payload  ratio,  normalized  deformed  shape  and 
stress  distributions  for  any  altitude  designated  as  BBAR  on  card  3.  Principal 
stresses  and  tension  field  patterns  are  computed  when  necessary.  A sample 
of  this  output  is  as  follows: 
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